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ABSTRACT 

Electrical prospecting is defined as the science and the art of determining the 
variations of the electrical constants (resistivity, magnetic permeability and the di- 
electric constant) of the earth’s crust and of interpreting these variations in terms of 
geological structure. The most successful systems of prospecting are based on the 
study of resistivity variations. The basic assumption made is that in general changes of 
resistivity follow the bedding planes. Electrical methods of exploration may be divided 
into two classes, direct current methods and alternating current methods. In part II 
the fundamental theory of direct current method is discussed and a typical survey is 
described. Part III deals with the theory of alternating current methods, with partic- 
ular reference to the optimum frequency to be used. It is shown that in general very 
low frequencies are desirable. Two alternating current surveys made by the Swedish 
American Prospecting Corporation are briefly described. In the conclusion (part IV) 
some of the difficulties of electrical prospecting are discussed. The depth to which in- 
vestigations may be carried is limited. In the present stage of the art, it would take 
exceptionally favorable conditions to obtain reliable information much in excess of 
2,000 feet. However, it is stated that improvements in methods of interpretation and 


in field technique should give electrical methods a definite field of usefulness in pros- 
pecting for oil. 


I. INTRODUCTION 


LECTRICAL prospecting may be defined as the science and the art of 

determining the variations of the electrical constants of the earth’s 
crust and of interpreting these variations in terms of geological structure or in 
terms of mineralization of a region. The three electrical constants whose vari- 
ation may be determined are the resistivity, the magnetic permeability and 
the dielectric constant. 

The resistivities of earth materials vary from a fraction of an ohm-cm 
to 10!’ ohms-cm. The materials having values below 100 ohms-cm are con- 
fined largely to mineral ores such as marcasite, galena, molybdenite, chal- 
copyrite and iron pyrite. In the study of oil structures the resistivities en- 
countered vary from 100 to 10'7 ohms-cm. In the Texas Gulf Coast region 
the common variation is from 200 to 20,000 ohms: cm with the exception of 
rock salt of which the resistivity has a value of the order of 10". 

The dielectric constant of earth materials has a value which ranges from 
1 to 10. The variation in magnetic permeability, except in rare cases, is 
negligible. 
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Systems of electrical prospecting are, therefore, confined to studying the 
manner in which the resistivity or the dielectric constant varies in the earth’s 
crust. The relative effect which these two constants have on the measurable 
electromagnetic quantities at the surface of the earth will be discussed later. 
However, it can be stated that up to the present time, the most successful 
systems of prospecting have been based upon a study of the resistivity varia- 
tions in the crust of the earth. 

The problems of any scheme of geophysical prospecting may be stated as 
follows: 

First, from mathematical and physical considerations it is decided what 
data must be taken at the surface of the earth in order to predict the dis- 
tributions of a chosen physical constant or constants in the earth’s crust. 

Second, instruments and field technique are developed to obtain these 
data. 

Third, mathematical tools are developed which translate the quantities 
measured at the surface of the earth into the distribution in the earth’s crust 
of the physical constant or constants under consideration. When these phys- 
ical and mathematical problems are solved there remains the problem of 
finding out what relation (if any) the distributions of the physical constants 
which are determined bear to geology and of interpreting the measurements 
in terms of geology. Furthermore, to be economically justifiable the overall 
cost from field measurements to translation into geology must be com- 
mensurate with the value of the information obtained. 

An ideal system of prospecting for oil would be one which directly in- 
dicated the presence of oil. Many attempts have been made to devise an 
electrical system which would meet this ideal requirement. However, these 
attempts have met with little or no success. In the present state of the art, 
the practical systems of prospecting for petroleum are those which are de- 
signed to give information concerning subsurface geology. 

The basic assumption in mapping geological structures by electrical 
methods is that in general changes of resistivity follow the bedding planes. 
Thus a clay bed or shale bed may have a lower resistivity than surrounding 
limestones or sandstones. These beds usually are parallel with or are only 
slightly inclined to the surface of the earth. In general, the lateral changes in 
resistivity are neglected, and the changes in resistivity with depth are deter- 
mined. If many such determinations are made over an area, the resistivity 
depth curves may be correlated and the subsurface topography may be ob- 
tained. 

Generally one or more beds which have a large resistivity contrast with 
the beds above or below them are used as marker beds in mapping geological 
structures. It is obvious that the bed which is being used as a marker must 
conform to the bed carrying the oil in order to be of any practical use as a 
marker bed in exploring for petroleum. In many localities the relatively shal- 
low beds are more or less parallel to the deeper ones. However, in many other 
places major unconformities occur and the shallow beds are no longer sub- 
stantially parallel to the deeper ones. In such cases it is necessary to follow a 
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marker bed which is below the unconformity and is conformable to the oil 
bearing horizon. The depth to which it is possible to work with any geo- 
physical scheme is, therefore, an important consideration. 

The depth which can be obtained reliably in electrical prospecting is a 
much debated question. Up to the present time the depth at which structures 
have been reliably mapped by electrical surveys has not exceeded 1500 feet. 
Electrical surveying has gone but a little deeper than the core drill. 

While electrical prospecting is thus a competitor of the core drill, it may 
be used as an aid in making correlations between core drill holes and also 
between wells in cases where geological correlation is difficult. This is ac- 
complished by making resistivity measurements in the holes at a series of 
depths and correlating on a resistivity basis. 

Electrical methods of geophysical exploration may be divided into two 
classes, direct current methods and alternating current methods. The first 
of these also includes the use of very low frequency alternating currents which 
are commonly used to eliminate the effects of polarization. The frequency, 
however, is low enough so that the effect of currents induced by the varying 
magnetic field may be neglected. Induced currents play the major role in the 
second, or alternating current method of prospecting. 


Il. Direct CurRRENT Metuops or PROSPECTING 


There are a great many direct current methods available for the determin- 
ation of the changes in resistivity with depth, but the same fundamental 


























Fig. 1. Lines of flow in a homogeneous earth. 


principles apply to all of them. A common practice is to pass the current into 
the earth by means of two electrodes on the surface, as shown in Fig. 1. If the 
earth is homogeneous, the lines of flow are arcs of circles connecting the elec- 
trodes. The potential on the surface of the earth is then called the normal 
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potential. Any changes of resistivity in the subsurface beds alter this normal 
potential. Thus measurements of potential on the surface of the earth give 
an indication of these changes of resistivity. Most prospecting methods differ 
only in the spacing of the electrodes and in the manner in which the potential 
is measured. The same mathematical theory thus applies to all of them. 

The theory is simplified by the fact that problems of steady flow in con- 
tinuous media are the same mathematically as similar electrostatic problems. 
A current J entering an earth of conductivity o at an electrode A may be re- 
placed by a charge of g=J/270 at A. The potential at any point in the earth, 
due to this electrode is V = J/2zer, where r is the distance from the electrode. 

If several electrodes A, are on the surface of the earth with currents [, 
entering the earth, the potential at any point P in the earth is 


1 
Vo— 2 (1) 
2ro Tn 
where r, is the distance from P to A,. It may be noted that algebraic values 
of the current must be used in this sum. Thus the potential at P in Fig. 1 
for the current flow between two electrodes is 


; I (- ~) (2) 
Dae ri ro j 


This additive property of the potentials is general and is not confined to a 
homogeneous earth. Thus we may simplify the theoretical discussion by con- 
sidering the potentials due to a single electrode. 

In order to illustrate the type of departure from normal potential that 
may be expected, a simple problem will be solved. The earth is of conductivity 
o,, to a depth a and is of conductivity a2, below this depth. A current J flows 
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Fig. 2. Coordinate system for two-layer problem. 


into the earth at the origin of cylindrical coordinates r, z, @. (See Fig. 2.) Let 
the potential in the surface layer be V; and the potential for z>a be V2. The 
conditions to be satisfied by the potential are: 
(1) The differential equations expressing the continuity of the currents 
av, 1 0Vi avy 
vV,=—_+- — = = 
Or? r or 03" 


VV. = 0. (3.2) 





0 (3.1) 


(2) The continuity of current and potential at the boundary 


V,=Vzatz=a (3.3) 
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o,0V,/dz = o20V2/dz. (3.4) 
(3) No current flow out of the surface of the earth 
oV;,/dz = 0 at z = O. (3.5) 
(4) The conditions near the electrode 
: I 
3-8 rg - 0.7 (3.6) 
2ra,(r? + 37)!/2 
(5) The conditions at ~ 
Vi-0as r->« 
(3.7) 


V.-O0O asryz-7 @, 


This problem has usually been solved by the method of images, which results 
in an infinite-series. An integral solution which is often more useful' may be 
obtained from the following solution of V?V=0 


V= i) [fe + fade Volar) dr (4) 
0 


where f; and fz may be any functions of \ subject to the condition that the 
integral converges. Thus if 


I «0 
Vv; = —— f [(A(A) — 1)e* + A(A)e™]Jo(Ar)dd (5) 
210; 0 
and 
V2= f B(d)e™*J o(Ar) dd (6) 
0 


where A (A) and B(A) are arbitrary functions of \ all the conditions (3) will be 
satisfied except the conditions at the boundary z=a. Of course, A(A) and 
B(X\) must be chosen in such a manner that the integrals converge. These 
functions are determined by substituting V; and V2 given by Eqs. (5) and 
(6) in the boundary conditions (3.3) and (3.4). The solution for A (A) is 


1 
A(\) = —————- (7 
) 1 + he-2h4 ) 
where 
02—~ Oj 
Be ee 
2+ 1 
On the surface of the earth, z =0, the potential is 
: I ) 1 — ke 2ra 
Vy, = — f —_——— | J,(ar)dd (8) 
210; 0 1 + ke she 


' This method is especially suitable when a number of horizontal layers are present. 
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I mk ae | 8) i (9) 
= — 0 ( 
2roir do LI + ke7e!r 
where 
B = »)r. 


When a/r is large, the integral becomes 


I ° 
— — Jo(8)dB = —— Ka (10) 


2mroir 2ro\r 


which is the normal potential for a homogeneous medium of conductivity 1. 


When a/r is small, the exponentials are approximately equal to unity, and the 
potential becomes 


I 1i-k I 
Pi nl ccm ete OE, (11) 


Irow 1+ k 2roor 


This is the normal potential for a homogeneous medium of conductivity o». 
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Fig. 3. Theoretical potential curves for two conducting layers. 


This example illustrates a fundamental property of this type of prospect- 
ing. The potentials at po nts further and further from the electrode are in- 








errr 





~ om & 











7 rt ee 








ELECTRICAL PROSPECTING 109 





fluenced more and more by the deeper layers. The surface layers become less 
important. In order to illustrate this effect values of u = 270,rV\/ J were plotted 
against values of r/a for various ratios of 02/0, (Fig. 3). The potentials are 
multiplied by 7 in order to eliminate the normal drop of potential with dis- 
tance. If the ground is homogeneous the potential curve is a horizontal line, 
the ordinate of which depends on the conductivity. Thus for an earth of con- 
ductivity a1, «=1 for all values of 7. For any other conductivity 02, « =0\/02 
for all r. In Fig. 3 the earth is of conductivity a, to a depth a, below which the 
conductivity is a2. The curves for u start at «=1 for r=0 then approach the 
lines u =0)/o2 asymptotically. These curves illustrate the gradually increasing 
effect of the lower medium as 7/a increases. 

Fig. 4 shows similar curves for a slightly more complicated case. There is 
a surface layer of conductivity o;, then a second layer of conductivity o2, and 
finally a medium whose conductivity is the same as that of the surface layer, 
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Fig. 4. Theoretical potential curves for three conducting layers. 


g,. The thickness of the surface layer is a and that of the second layer 2a. Two 
curves are shown in Fig. 4: (a) for ¢.=30,, and (b) for a; = 302. In curve (a) 
u=1 at r=0, and then wu decreases for increasing r reaching a minimum at 
about r= 2.5a where the second layer has most effect and finally approaching 
the line «=1 asymptotically. Curve (b) rises from «= 1 when r=0 to a max- 
imum and then decreases rapidly as it approaches its asymptote u=1. Thus 
the effect of the second layer on the potential is small near the electrode. It is 
a Maximum at a distance equal to roughly the depth of the layer, and the 
effect diminishes as r increased beyond this point. These simple cases are 
sufficient to show how the potential due to a single electrode is affected by 
vertical changes in conductivity. 

In practice it has usually been found more convenient to measure po- 
tential differences rather than absolute potential. One method is to put the 
current electrodes a long distance apart and measure potential differences in 
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the neighborhood of one electrode so that the potentials due to the distant 
electrode may be neglected. The potential difference is measured between two 
search electrodes at distances 7 and 2r from the current electrode. If the earth 
is homogeneous and of conductivity ¢, the potential difference between the 
search electrodes is: 


Iri 1 I 
Ve = Vy) — VU) = —|- = =| oS (12) 


2roLr 2r 4dror 
The resistivity p may, therefore, be determined as 
p = 1/o = 4nrV,/I. (13) 


If the earth is not homogeneous, Eq. (13) may be used to define the apparent 
resistivity. Let V,(r) be the potential of the current electrode. Then the ap- 
parent resistivity is 


4rr 
a 2 [Vi(r) — V,(2r)]. (14) 


This function is not independent of 7 except in the special case of a homo- 
geneous earth. The variation of the apparent resistivity with r is used to de- 
termine the changes of resistivity with depth. As may be expected from the 
potential variation, the apparent resistance is equal to the resistance of the 
upper layer when ¢ is small. As 7 increases the apparent resistance approaches 
the resistivity of the second layer. Long before it reaches this value, a third 
layer may affect it, and then the apparent resistivity will have a tendency to 
approach the resistivity of this layer. When a great many layers are present 
the resultant curve may be rather complicated, and the interpretation diff- 
cult. 

A more common method of prospecting which leads to the same value for 
the apparent resistivity was developed by Wenner and Gish and Rooney. 
The two current electrodes and the two search electrodes are placed in a line 
equidistant from one another. 


‘—foe—reoeere-: 


A Cc D B 


The current enters the earth at A and leaves at B. The potential difference 
is measured between C and D. The common spacing is r. The apparent re- 
sistivity is again given by Eq. (14), where V; again represents the potential 
due to a single electrode. The advantages of the method are that the electrode 
spacing need not be excessively great, and that possible errors due to neglect- 
ing part of the field do not enter. However, to obtain values of the apparent 
resistivity for different values of 7, all four electrodes must be moved instead 
of only two as in the previous method. This method is often called the Gish- 
Rooney method of electrical prospecting. 

In order to illustrate the results that may be expected under favorable 
conditions, a typical Gish-Rooney survey which was taken along a line of 
core drill holes is shown in Figs. 5 and 6. Fig. 5 shows the curves of apparent 











: 














Ee 























ELECTRICAL PROSPECTING 111 





resistivity at the six stations shown on the profile of Fig. 6. These curves 
show a layer of very high resistivity which is at the surface at stations A and 
B on the west, and then dips to the east. A layer of low resistivity is on the 
surface at stations C, D, E and F, below which there is the high resistivity 
layer. There is apparently a very thick layer of low resistivity below this high 
resistivity layer. A rough qualitative interpretation of these curves is shown 
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Fig. 5. Apparent resistivity curves obtained by field measurements along profile of Fig. 6. 


in Fig. 6. It may be seen that it checks with the core drill holes which show a 
monocline dipping to the east. 

Another method of utilizing the apparent resistivity, the ground resistiv- 
ity map, has been developed by Schlumberger. An electrode spacing is chosen 
which is approximately equal to the desired depth of investigation. The ap- 
parent resistivities with this electrode spacing are measured at various points 
on the surface of the earth. A contour map is then made of these values. This 
method is particularly suitable in a country with a shallow overburden which 
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covers a thick sediment with a large change in resistivity. Then the contours 
will roughly approximate the topographical contours on the sediment. Sup- 
pose that the sediment has a low resistivity. Then when it is near the surface, 
the apparent resistivity will be low. When the sediment is deeper, the ap- 
parent resistivity will be higher. 

The apparent resistivity map, of course, gives only a qualitative inter- 
pretation of the electrical data. The methods at present in use for interpret- 
ing Gish-Rooney surveys give only approximate depths to the points at 
which sharp changes occur in the resistivity of the earth. These approxima- 
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Fig. 6. Qualitative interpretation of Gish-Rooney survey shown in 
Fig. 5, giving correlation with core drill. 


tions become worse and worse as the number of resistivity contrasts become 
greater. A rigorous solution of the problem of finding the variation of re- 
sistivity with depth from surface potential measurements under the assump- 
tion that the resistivity is substantially a function of the depth only could be 
based upon solutions of the potential problem given by the theory of images. 
However, these solutions are cumbersome to handle when the number of 
layers which must be dealt with exceeds two. The integral solution given in 
this paper is much easier to handle in the case of multilayer problems and 
should aid in placing the interpretation of direct current surveys on a more 
exact basis. 


III. ELEcTROMAGNETIC METHODS 


Electromagnetic methods depend primarily on the induction of currents 
in the earth by means of a primary current on the surface. Inhomogeneities 
in the earth will affect the resultant field on the surface. In oil prospecting 
these inhomogeneities take the form of practically horizontal beds with 
different conductivities. Many different schemes have been used with vary- 
ing success. The factors which should be taken into consideration are: 

(1) The primary field should be as simple as possible so that the associated 
theoretical problems may be solved. 

(2) The frequency should be chosen so that the inhomogeneities which are 
being followed will have a maximum measurable effect on the field. 

(3) The field components should be measured at points where the beds 
produce the greatest measurable anomaly in the field components. Thus 
both the magnitude of the change and the percentage change in the field 
should be considered. 

The quantities which may be measured are the magnitude and time phase 
of the components of the electric and magnetic vectors on the surface of the 
earth. 
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There are two simple primary fields for which the theoretical problems 
have been solved; the field of a circular coil, which may be considered to be an 
oscillating magnetic dipole with its axis perpendicular to the plane of the coil 
and the field of a long horizontal wire. Most electromagnetic methods have 
made use of one of these. 
The basic equations that govern alternating currents are Maxwell's 
equations. The field components must satisfy the wave equation 2 
; a4, A, ; 
V?7A, — Ss orien = 0). (15) 
ot ot? 
If the impressed forces are simple harmonic functions of time, A; can be 
represented as 


A, = real part of Ae~™! 
where A is in general complex. Then Eq. (15) becomes 
V7A + iwpA(o — wp) = 0. (16) 


This equation shows the relative importance of the conductivity and the rel- 
ative permittivity or dielectric constant in determining the current flow in 
the earth. If c>wp the conductivity will be the determining factor. If wp>a 
the dielectric constant will be the determining factor. If the conductivity is 
10-‘ohms-cm and the dielectric constant is 10, then the two are of equal im- 
portance when the frequency is 


f = ——__—___—_——_ = 1.8 x 10’. 
2x X 10 XK 8.85 XK 10-4 

It will be shown later, however, that the high frequency waves are very highly 

damped by the surface layers, and are thus not important in oil prospecting. 

At the lower frequencies, the conductivity governs the distribution of current 

in the earth. Methods of locating oil directly by means of its dielectric con- 

stant are thus not very hopeful. 

In order to show the importance of using a low frequency, a particularly 
simple situation has been assumed. The primary field is due to a vertical oscil- 
lating electric dipole on the surface of the earth.’ The earth is homogeneous 
and of conductivity ¢ to a depth h, below which there is a layer of very high 
conductivity which we shall assume to be infinite. The problem to be solved 
is—at what frequency will the highly conducting layer have the biggest in- 
fluence on the field at the surface of the earth? 

The field may be conveniently represented by the Hertzian function I, 
which satisfies the wave equation 

2 The rationalized practical system of units is used. Magnetic intensity H/ is expressed in 
ampere turns per cm. Magnetic flux density B is expressed in webers per cm?(=10* Maxwell's 
per cm?*). Permeability u=B/H=42 X10~° for free space. Electric intensity E is expressed in 
volts percm. Displacement D isexpressed in coulombs per cm?. Permittivity p= D/E =8.85 +107" 
for free space. Conductivity ¢ is expressed in mhos per cm. 

? Riemann-Webers, Diff. Gleich. der Phys. (Seventh ed. 1925) p. 542. 
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; : jh = ko? = wpow*® = w*/c? in the air , 
Vl + kil = 0 - (17) 
lk = k= pw? + ivow in the earth 


| 
' 


where c is the velocity of light. 
Cylindrical coordinates (7, z, @) with the origin at the dipole will be used. 
The electric intensity E and the magnetic intensity H/ are: 


E; = H,= H, =0 
E, = 0*Il/drdz 
; 1 @ oll 
E,=-— — r=) (18) 
r Or Or 
ik,? dll 
Hw=—-—- 
ww Or 


We shall take the Hertzian function for the air to be IIp, and the Hertzian 
function for the earth to be II. The boundary conditions then are 
ko*IIo = ker | 
e - > atz=0 (19) 
All,/dz = dx/dz § 
and 
oll 
—=Oatz=—Ah. (20) 
02 


The Hertzian function for a dipole in free space is 





x Xr : — 
Hy = | 20 hu?)"" F(A) dd, (21) 
0 (2 = ko*)! > 


We, therefore, take for I, 


" d xiii 
IIo = f | = —— fo) eke)" Jo(Ar)dd 2z>0 (22) 
0 (\? _ k*)! 2 


and for II 








I = f [gr(AjezO—F)'* 4 go(Nem2*- FT (Ar) aX. (23) 
0 


The functions f, gi, and gs are chosen so as to satisfy the boundary conditions 
Eqs. (19) and (20). The solution for Ip is 


* 2k2r|W2 Oho?) TOA) dd 
IIo = f . : , (24) 
o RAZ — Ro2)2 + ko2(A2 — &2)!/2 tanh A(A2 — &2)2/2 





The normal field for a homogeneous earth of conductivity ¢ is 


I, = 





f 2h2reW FO) TOA) ad 
0 


k2(N2 — ko?)/2 + bo2(A2 — &2)1/2 , 
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The change in the field due to the layer of high conductivity is: 
ly — I, = f eapero 1 soanyan | (26) 
0 
where 
2k*ko?(A2— k*)(1—tanh h(d?— k? tail *) 


eeensaisnaiiaa ———— (27) 


)= [R2(A2— eo?) "2 ho2(A2— R2)"/2] [ h2(A2— Ro?) 2 ho A2— 2) "tanh h(A2— k2)"2] 








g(r 


In order to calculate the effect of the layer on the normal field, one could cal- 
culate the components of this change in field and compare them with the 
normal field components. A simpler way to compare them is to compare the 
power flow associated with the change in field to the normal power flow. The 
changes in the field components are: 


(E,)o — (Ey)n = J P(A)(A2 — ko?) P/2e-2O*— ho) ZT (Ar) Ar 
ik (28) 
tRo 

(Hs)o — (Hs) = = (rez KY" NET AAP)AD. 
pw 

The power flow of this field through a horizontal plane in the air is: 
P=r J [(Es)o — (E,)n] [Heo — (He)a]*rdr 
(30) 


ark? 





— (ko? — d2)"/2| H(A) | 2A5dd. 


The asterisk indicates that the complex conjugate of (J74)y9—(/7s), is to be 
taken. 

For all but very high frequencies, k? is very much greater than k,o?. Under 
this condition, the power flow is approximately 


4rko*[(1 — R)? + 1°] 1 


" Nit saitamanes 31 
(2)"%yolkI(1-R+D °° R-T (1) 





where R and 7 are real, and 


i+: 


J1/2 





tanh h | =R+il. (32) 


The normal power flow from a homogeneous earth of conductivity @ is 





—__—_—__—_ (33) 


4rko?| k| 4 pt (ko? — d?)*/2\3dd 
7 ——— f 
o | BAZ — be)? + y2(A2 — B2)12| 2 


pw 
which becomes for k?>k,? 
P, = Srko’/3p. (34) 
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The ratio of the power flow of the change in field to the normal power flow is: 


p/P 3 Ro 2.81 XK 107? , 
_=— ——— gu) ————_ g (8 (35 
2 22) k\ w sali ho s\#) } 
where 
(1 — R)? + I? 1 
g(s) = « ———— log ———- (36) 
i—-R+] R-I 
u = h(pow)'? (37) 
1+ 
tanh «w——- = R+il. (38) 



























































Fig. 7. Power ratio curve for vertical electric dipole. 


Fig. 7 gives a graph of g(u) which shows a maximum at u=0.29. This gives 
the optimum frequency in terms of the depth of investigation and the con- 
ductivity of the upper beds. This frequency is 

10°u? = (0.29)?-10® 10? 


—— = (39) 
Saha h,,*6 


S27 ho 
if h is measured in meters. If ¢=10~* and 4=100 meters, the best frequency 
to use is about 100 cycles. If h = 300 meters the optimum is about 10 cycles. 
At higher frequencies the effect of the conducting layer on the field becomes 
smaller. Thus for ¢=10~* and h=100 meters, the effect of the layer with a 
frequency of 2500 cycles is only about one tenth of the effect at 100 cycles. 
At very low frequencies only a small percentage of the power input enters 
the ground. At high frequencies, most of the power that enters the ground is 
absorbed by the surface layers. 

A similar investigation has been carried out to determine the effect of a 
highly conducting layer on the field of a long horizontal wire carrying an 
alternating current. As in the previous case, the earth is assumed homogene- 
ous with conductivity o to a depth h, below which the conductivity is infinite. 






















ELECTRICAL PROSPECTING 117 
























The wire is on the z-axis of rectangular coordinates, the y-axis is vertical 
and the x-axis is horizontal and perpendicular to the wire. Maxwell's equa- 
tions for the conducting medium reduce to: 


oll, oH, 


— =(o — wp)E 
Ox Oy 
(40 
0k, ; 
= twull, 
Oy 
Ok, : 
— =—_—— wu , 
Ox 


V7E. + k*E, = 0. 


The equations for the air may be obtained by setting ¢=0 and p= po. Let us 













: take 
E, = Eo in the air 
iy (41) 
FE, = E in the earth. 
The boundary conditions then are 
Ex, = E ) 
Ok dV = OE /dy 
kE=Oatyvy=—h (43) 
and 
Ey — ko log (x? + y?)'/? as x, vy 0. (44) 
The solution of these equations for Ey is 
. " 2e~¥™ cos Axdd 
Eo = ko f eee 4 cpa a ‘ ( $5) 
o (A? — ko?)!/2 + (A? — &*)!/? coth h(? — k?)'/? 
The normal field for a homogeneous earth is 
' : oeiil 
_ ° 2e~ vi —k)"" cos Axdd 
E, = ko f — (46) 
Jo (A — ko?)!/? + (A? — R?)1/? 
Thus the change in the field due to the conducting layer is 
Ey, — E, = ko [ o(r)em¥ =o") "* Cos Kxdd (47) 
/7f 
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where 


‘ 2(d? — k?)*/2(1 — coth h(d? — k?)*/2) 
(A) = wah 7 a 
[ (A? — bo?)!/2-+ (A? — h?)"/?] [(A? — ho?) /2-+ (A? — ?)!/2coth h(A? — k?)"/?] 








(48) 


and the change in magnetic intensity is: 
iko (* » ii tire 
(H.)o — (Hz), = — (A? — ko)*o(A)e~4¥® —*” "cos Axdd. (49) 


MW J/0 


The power flow through a horizontal plane in the air per centimeter of wire 
length is 


P=} [Eo — Ex) [(H ~ en a "| g(a) | Ube? — 8)"2aa. (50 
2 “90 — En z)0 ne P(A) | “(Ro )©"dKX. (99) 


When k*>&,* this integral reduces to 





r*kot 
a — |[(1 — R)? + /?| (351) 
duw| Rk}? 
where R and 7 are real and 
i+i 
R+ i = tanhh F\ ( a): (52) 


The normal power flow through the plane for a homogeneous earth is: 


1 ko! i 
> = —___.. (53) 
duw | Rk?) 


As the frequency approaches zero, the ratio ko'/w k |? approaches zero, so 
that both the normal power flow and the power flow of the change in field 
due to the layer approach zero. The ratio, however, is 

P/P, =(1—R)?+7? (54) 
where 
(1 + i) 
(2)"2 
u=h\ k| = 0.0281 Munreter(of)!/? 


R+ il = tanh u 





The ratio p/p, is shown graphically in Fig. 8. It is a maximum at zero fre- 
quency, and then drops off as the frequency increases, approaching zero at 
very high frequencies. In this case the optimum frequency depends primarily 
on the measuring apparatus. The frequency should be chosen as low as the 
apparatus permits the field quantities to be measured accurately. 

The curve in Fig. (8) shows the absorption of the surface layer. If the 
conductivity of the surface layer is 10-* ohms and the depth to the conducting 
layer is 100 meters, the effect of the layer at a frequency of 3000 cycles is only 
























ELECTRICAL PROSPECTING 





119 


about one twentieth of the effect at 100 cycles. This curve illustrates the im- 
portance of using a low frequency alternating current. 
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Fig. 9. Electromagnetic survey across a salt dome. 
In order to illustrate the results that may be obtained in the field by a 


good electrical survey under favorable conditions, two surveys made by the 
Swedish American Prospecting Corporation for the Gulf Oil Corporation will 
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be described. The method of survey used was the Sundberg method,‘ which 


will be briefly described. 


The primary field is the field of a long horizontal wire carrying an alternat- 
ing current with a frequency of 100 to 500 cycles. Measurements of the direc- 
tion and time phase of the magnetic intensity vector are made at a series of 
points along transverse lines. These measurements are made by means of a 


search coil connected to a compensating arrangement. 
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Fig. 10. Fault traced by electrical methods where only one well was drilled in the field. 


The induced currents are assumed to be confined to a series of beds which 
have a much higher conductivity than the surrounding beds. Each of these 
beds is characterized by a certain induction factor p which depends on its 
conductivity and thickness, and also on the frequency of the exciting current. 
From the measurements made at the surface, the induction factors and the 


* T. Zuschlag, Tech. Pub. No. 313 Inst. of Min. and Met. Eng. 
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depth of the beds are determined. The conducting beds which are followed 
are called marker beds. 

Fig. 9° shows the results of a survey across a shallow salt dome. Three con- 
ducting beds were assumed, one at the surface, one at an intermediate depth, 
and the last at about 400 feet. The intermediate bed is arched over the dome. 
The deeper conductor is missing where it is pinched out by the salt, and 
locates the flanks of the dome. 

Fig. 10 gives the electrical indication of the Bruner fault zone, Texas. A 
clay bed at about 500 feet was used as a marker. When the survey was started 
there was only one producing well in the area. Later production has followed 
parallel to the electrically indicated fault. The electrical picture was not only 
qualitively right, but later drilling showed a constant interval between the 
clay bed and the Austin chalk, which is the producing horizon. Thus, a very 
good shallow picture of this fault was obtained electrically. 

IV. ConcLusIoNn 

Up to the present point we have discussed in a general way the problems 
encountered in electrical prospecting under very simple geological conditions. 
Experience in the field has shown that in general good marker beds can be 
found within the depth to which investigation can be carried electrically. 
However, in countries where the shallow sediments have been laid down in 
the form of lenses, these beds may not be continuous enough to give a struc- 
ture picture. 

On monoclinal slopes where successive beds are outcropping so that the 
surface material is changing in character, the lateral variation in conductivity 
of the surface material is large. This resistance variation of the near surface 
material makes the problem of interpreting electrical measurements some- 
what difficult. A condition of this kind was encountered in certain parts of the 
Balcones fault zone of Texas. The variation in surface resistivity resulted in 
many apparent fault indications and it was difficult to distinguish these ap- 
parent indications from the true ones. 

The depth to which geophysical investigations can be carried is of con- 
siderable importance in regions where the lower beds which contain the oil 
are not laid down conformably with the shallow beds. If the beds in which oil 
must be sought lie below an uncomformity which occurs at a depth of 5000 
feet, then it is obvious that a prospecting scheme cannot be used which does 
not give information, below 5000 feet. 

Many extravagant claims have been made for the possible depths to 
which it is possible to work by electrical methods. In direct current cases, if 
the earth were homogeneous except for two or three sharp breaks in conduc- 
tivity, it would be possible to map to a large depth. However, even in this 
ideal case as the depth of investigation is increased, there would be a gradual 
loss in detail. Actually conditions are much more complicated. As the depth 
of investigation is increased, more and more discontinuities enter both later- 


° Surveys shown by Figs. 9 and 10 have previously been published by Zuschlag in 
A.I.M.M.E. Technical Publication No. 313. 
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ally and vertically. Consequently, the difficulties of interpreting the measure- 
ments become more and more difficult. Up to the present time information of 
immediate value in oil geology has not been obtained below a depth of about 
2000 feet. 

In the alternating current case there is a definite depth of current penetra- 
tion for a given frequency and a given earth resistivity. In attempting to work 
at deeper and deeper depths, the operating frequency must be lowered con- 
tinually. In addition to the difficulties mentioned in the direct current case, 
one is confronted with the problem of measuring very weak electric and 
magnetic fields of low frequency. In some cases where the resistivity of the 
surface materials is high the alternating current field is less influenced by 
lateral irregularities in the surface layer than the direct current field. In the 
present state of the art, it would take exceptionally favorable conditions to 
obtain reliable information at a depth much in excess of 1500 feet. 

At the present time electrical methods of prospecting for oil seem to be in 
disrepute. This is partly due to cost of electrical surveys as compared with 
other geophysical methods and partly due to the failure of the extravagant 
claims made for the process to materialize. However, the electrical method of 
prospecting for oil cannot be forgotten because it is one of two prominent 
geophysical methods in which it is possible to control the field being em- 
ployed. Improvements in methods of interpretation and in field technique 
should give electrical methods a definite field of usefulness in prospecting for 


oil. 
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A NEW INSTRUMENT FOR MEASURING VERY 
SMALL DIFFERENCES IN GRAVITY 
By KenNetu HartLey 
Hartvey Gravity BALANcE Corporation, Houston, Texas 
(Received January 2, 1932) 


ABSTRACT 

Description of a new portable instrument for measuring relative values of gravity 
to within two or three parts in ten million, specially designed for geophysical explora 
tion. The principles of the design are analysed and the methods for eliminating effects 
of elastic hysteresis, temperature changes, variations in barometric pressure, etc., are 
discussed. Also effects of initial stresses in materials, defects in alignment of locking 
mechanism, inaccurate leveling, etc. These difficulties are serious but seem to have 
been overcome. Results of preliminary field measurements near Houston, Texas, are 
presented and some comparison between the type of information given by this instru 
ment and by the torsion balance. 


NTIL now the only successful method of measuring gravity has been 

with some form of pendulum, in other words determining the accelera- 
tion produced, or, we may say, obtaining the value of gravity in terms of 
inertia. There is only one other possible method and that is to measure it in 
terms of the elasticity of a spring which can be calibrated in some way. This 
has the advantage of eliminating time measurements but introduces some 
other difficulties that are avoided by the pendulum. 

In all previous attempts that I know of to use a spring the effort has been 
to obtain extreme sensitivity by some arrangement approaching unstable 
equilibrium, and one reason for their failure is that this unstable condition 
also exaggerates all the causes of error and uncertainty in the indications of 
the apparatus. In the new instrument the effort has been in the other direc- 
tion, that is to obtain the maximum stability, the small displacements are 
then amplified by optical means with less uncertainty. 

The fundamental idea of the design is that if we can apply to a suspended 
mass an upward force almost, but not quite, equal to the pull of gravity and 
maintain this force constant then it will be easy to measure changes in the 
small additional force required to hold the system in equilibrium. The first 
problem then is to apply a force that can be maintained constant to the re- 
quired degree of accuracy. The value of such an instrument for geophysical 
prospecting depends on its ability to measure very minute differences in the 
gravitational force, so the instrument has been designed to have the greatest 
sensitivity that can be of any practical use. Since the change in the gravita- 
tional attraction due to a difference of one foot in elevation is almost exactly 
one part in ten million this was taken as the maximum sensitivity that could 
possibly be useful; for most purposes three or four parts in ten million will be 
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all that is desired. This is about ten times the accuracy of the best work done 
with the pendulum apparatus of the United States Geodetic Survey. 

Figure 1 is a simplified diagram of the essential parts of the apparatus. The 
main spring, which carries about 99.9 percent of the load, is formed from 
heavy wire and not stressed to more than 20 percent of its apparent elastic 
limit, so as to give the greatest possible dependability. A very light beam, 
hinged at the left side, is flexibly attached between the spring and the weight 
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Fig. 1. Schematic diagram of gravimeter. 


and carries two small mirrors at the free end. These mirrors act on the prin- 
ciple of the optical lever and are so mounted that a vertical displacement of 
the beam causes them to turn in opposite directions. This doubles the ampli- 
fication but a much more important fact is that it avoids the necessity for 
maintaining a fixed point of reference in the telescope; when the two mirrors 
are in the same plane the two images of the lamp filament as seen in the eve- 
piece will be in one straight line no matter how much the eye-piece and other 
optical parts may be displaced. A very light weighing spring is attached to the 
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beam between the center and the hinged end and the tension in this spring 
is controlled by a micrometer screw above. 

If the instrument is adjusted at a certain location so that the two images 
of the lamp filament are in line and it is then taken to a different location, the 
difference in gravity will cause a slight displacement of the images which are 
then brought back to the zero position by altering the tension in the weighing 
spring; the difference between the two readings of the micrometer dial is 
then a measure of the difference in the force of gravity at the two places. 

The difficulties involved in carrying out this idea to the necessary degree 
of precision have been greatly increased by the requirement that the instru- 
ment be rugged enough to stand transportation over rough roads and small 
and light enough to be carried by two men where there are no roads. 

Suppose a mass of 100 grams or more suspended by a spring of such stiff- 
ness that the total elongation, from no load to full load, is 100 millimeters; 
then, to maintain the tension constant to one part in ten million would re- 
quire; first, that the length be controlled to within one one-hundred-thou- 
sandth of a millimeter; second, that the temperature be maintained constant 
for the whole period required for a series of observations so that the elonga- 
tion of the spring due to temperature effect minus the linear expansion of the 
supporting frame will be within the same limit; third, that effects of elastic 
hysteresis be eliminated not only from the main spring but from the flexible 
supports of the beam and mirrors also. 

In the present instrument the distance between the thin ribbons sup- 
porting the mirrors is one millimeter and the distance from the mirrors to the 
eye-piece is 600 millimeters, so that the amplification for each mirror is 1200 
times, but as the two mirrors move in opposite directions the relative dis- 
placement of the images is amplified 2400 times. As the beam ratio is 1.8 the 
total amplification is 4320, but as this is observed through an eye-piece with 
a magnification of 14 the apparent amplification is over 60,000, so that a 
movement of.the suspended mass of one one-hundred-thousandth of a milli- 
meter will look like six-tenths of a millimeter and be easily observed. If the 
force of gravity is increased by one part in ten million, part of the increase is 
used in stretching the main spring and part in bending the small filaments 

that support the beam and mirrors, so the sensitivity in terms of force is not 
so great as this, but the first requirement is accomplished—the maintainance 
of the length of the main spring to within one part in ten million. 

I expected to be able to make the beam and mirror system so flexible that 
the resistance to movement would not be any greater than that of the main 
spring, but I have not succeeded in doing so. Although the filaments that carry 
the mirrors are only 0.01 mm thick, it is the stiffness of these filaments that 
limits the sensitivity of the instrument. The present instrument is satis- 
factory for readings of one part per million but is difficult to read any closer 
than that so a new beam system is being made which is expected to be at 
least four times as sensitive as the present one. 

The temperature effect on the spring is very large and unless accurately 
controlled will completely mask the smal! variations in force that are to be 
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measured. The ideal material for the spring would be a nonmagnetic metal 
with a high elastic limit but a low value for Young’s modulus and a low tem- 
perature coefficient. This combination might be found in one of the less com- 
mon metals or alloys but information on the elastic properties of these 
metals, particularly temperature effects, was not to be had. Assuming that 
the temperature coefficient of elasticity would probably be smaller in a metal 
with very high melting point led to the consideration of tantalum, molyb- 
denum, etc., and several springs were made of approximately the dimensions 
that would be required and tested in that form. The final choice was an alloy 
of tantalum and tungsten which seems to be fairly satisfactory. The spring is 
wound with a high initial tension so that when carrying its full load the coils 
are separated only about half a millimeter, making the total length of the 
stretched spring less than ten centimeters. The frame of the instrument is 
made of an aluminum alloy, which has the highest coefficient of expansion of 
any suitable metal, and is constructed on the principle of the old grid-iron 
compensated pendulum so that the total length of the high expansion 
members is about 80 centimeters, the low expansion rods are invar steel. This 
compensation is mainly for the purpose of avoiding too great a change in the 
relative position of the parts when the working temperature is altered. In 
order that the whole instrument be in a condition of thermal equilibrium and 
free from temperature stresses it must be kept at a constant temperature day 
and night during a whole series of observations; it is therefore enclosed ina 
new type of portable thermostat, designed on the principle worked out at 
the Bureau of Standards two or three years ago, which controls the tempera- 
ture to within one ten-thousandth of a degree. The heating is electrical and 
is designed to work with a six volt storage battery. 

A study of the action of this instrument while the temperature is either 
rising or falling will quickly convince any one that any plan for making a 
temperature correction, instead of maintaining a constant temperature, is 
wholly out of the question. In the first place, even with the best possible com- 
pensation the temperature would have to be measured to within one hun- 
dredth of a degree, and in the second place, it is impossible to have the entire 
apparatus at the same temperature unless it is kept constant for a consider- 
able time. In the present instrument if the temperature is altered by two or 
three degrees it requires at least five hours to reach a condition of equilibrium 
and if the change is much greater it may take all day, but with thermostat 
control the temperature effects are also brought within one part in ten million. 

The change in density of the air, due to varying barometric pressure, may 
alter the buoyant force on the suspended mass by several times the quantity 
that we are attempting to measure, so the case is closed air-tight and the two 
operating rods are surrounded by a new type of mercury seal which prevents 
the passage of any air in or out even with a difference of pressure of more than 
two inches of mercury. 

The effects of elastic hysteresis are eliminated by locking the weight and 
spring in the zero position between observations, so that the tension on the 
spring is never varied by more than the difference between the force of gravity 
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at two successive stations. Any hysteresis effect is then only a small fraction 
of this amount. This sounds simple but it is difficult to carry out for it re- 
quires that the suspended mass be locked to within less than one thousandth 
of a millimeter of the true position. To meet this requirement and to avoid 
any strain on any part of the spring from jolts in transportation has neces- 
sitated a very elaborated locking mechanism. 

In connection with this is another matter that no one seems to think of 
unless he has been working with a seismograph; that is vibration. The ground 
is constantly vibrating, and a weight hung on a spring is an excellent seismo- 
graph if its motions can be closely observed, so some way must be found to 
reduce the sensitivity to vibrations while maintaining the sensitivity to 
changes in the static force. In the present instrument this is accomplished to 
a certain extent by excessive air damping but this is only satisfactory under 
very favorable conditions so some modifications are being planned for the 
beam system which will give a more favorable ratio between these two effects. 
In the city traffic vibration is the chief difficulty and in the country it is the 
wind in the trees, so it is impossible to use the instrument in very windy 
weather but with a light wind readings can be made by watching for a favor- 
able moment, the mirrors may be almost stationary for several seconds at a 
time and an experienced operator can utilize these moments to secure read- 
ings to about one part in a million. 

One more difficulty encountered in the construction of the instrument 
was due to initial strains in the frame. Any change in temperature would 
alter the elastic forces present and cause the frame to twist enough to throw 
the beam system badly out of adjustment. To avoid this as far as possible the 
beam and mirror system is carried by a separate frame entirely independent 
of that which carries the rest of the mechanism. This frame is made of an- 
nealed aluminum tubes very accurately fitted into sockets in the top and 
bottom plates, and clamped in position by a method designed to avoid put- 
ting any strain in the tubes when they are tightened up. 

It was found, in the first experimental instrument, that the temperature 
stresses together with stresses due to imperfect construction of the locking 
mechanism would pull the beam and mirrors slightly out of alignment when 
it was locked. This caused distortion of the filament supports of the mirrors 
and of a light guide spring that was used in that instrument and when re- 
leased the filaments would not immediately return to their relaxed condition 
but there would be a slow drift of the readings extending over a period of 
from half an hour up to several hours before the correct reading was obtained. 
In the latest instrument this effect is almost eliminated but it is sometimes 
necessary to wait about ten minutes. 

These effects due to elastic hysteresis have been the most serious diffi- 
culties encountered, and there is still room for improvement in this direction, 
but measurements are now being made than can be checked, day after day, 
to +1-10-* and, as soon as a small alteration in the mirror system can be 
completed, which will give a better ratio between the effect of gravity and the 
effect of vibrations, I expect to get an accuracy four or five times as great. 
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The instrument must be very accurately leveled, that is, the axis of the 
instrument must be in the direction of the gravitational force. An angle of 
ninety seconds between the axis of the instrument and the true vertical would 





Fig. 2. The gravimeter controls and adjust ments. 
theoretically cause an error of one part in ten million, but it would cause a 
much greater difference in the actual measurement because of lateral strains 
in the filaments carrying the beam and mirrors. The present instrument is 





Fig. 3. The gravimeter in the field. 


easily adjusted to thirty seconds and that seems to be sufficiently close for 
the purpose. 

The dial readings give a direct measure of the differences in gravity and 
do not require any correction except to multiply by the calibration factor for 
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the instrument and subtract the effect of elevation, which is one part in ten 
million for every foot. These differences in elevation can be determined to 
within three or four feet with the Paulin Altimeter, provided that it is used 
according to a particular system that has been worked out for this purpose, 
but if greater accuracy is desired the land must be surveyed. These facts 
make the use of the gravity balance much simpler than any other geophysical 
instrument except the magnetometer. 

Recent field work has shown that in reasonably good weather it is easy 
to make three stations per hour. The torsion balance makes three stations per 
day and the Government pendulum apparatus three or four stations per 
month. This measurement which takes only a few minutes appears to have 
four or five times the accuracy of a single run of the pendulum which takes 
twelve hours. The information obtained is exactly the same as with the pen- 
dulum, the total intensity of the gravitational field. For this reason the pres- 
ence of a heavy mass near by on the surface or at small depth has no effect on 
the reading; the presence of such a mass merely alters the direction of the 
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Fig. 4. Gradient and gravity curves for assumed structure. 


resultant force and as the instrument is set exactly in line with this direction 
the measurement is the same. These surface masses make so much trouble 
for the torsion balance operator that it seems to be hard for him to believe 
that they do not disturb this instrument. 

The relation between the gravity anomaly, which is measured directly 
by the gravity balance, and the gravity gradient, which is measured by the 
torsion balance, is very clearly shown in a diagram in Ambronn’'s Elements of 
Geophysics, which he has taken from a paper by von Eétvos. Figure 4 is 
derived from this diagram; the heavy curve representing the values of the 
gravity gradient and the light curve the gravity anomaly caused by a mass 
of heavy rock at varying depths below the surface. It will be noted that the 
curve representing the magnitude of the anomaly follows very closely: the 
outline of the rock mass below while the curve representing the magnitude of 
the gradient bears no resemblance whatever to this outline and requires some 
study to interpret its meaning. In actual field work these curves are not ob- 
tained but only a rather small number of points on them. A comparatively 
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small number of gravity measurements would give the anomaly curve with 
a fair degree of accuracy but the same number of ordinates to the gradient 
curve would leave it still extremely uncertain and the results consequently 
very difficult to interpret and the interpretation of very doubtful validity. 

While the new instrument is still in the experimental stage and there will 
undoubtedly be improvements in details of construction, it is not probable 
that it can ever be used effectively by an unskilled operator. The method of 
operation is extremely simple but it seems to be necessary for the operator 
to have had some kind of experience with precision measurement in order to 
guard against deceptive appearances and be sure that the instrument is 
really in equilibrium when the reading is made. 
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CHARTS FOR TORSION BALANCE READINGS* 


By M. M. Stornick 


Geopuysics DEPARTMENT, HUMBLE OIL AND REFINING COMPANY 
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HE calculations necessary to determine the values of the gradient and 
curvature quantities from the readings taken at a torsion balance station 
involve either slide-rule or logarithmic manipulations which are somewhat 
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tedious, if not entirely cumbersome. It is quite a simple matter to show that 
the result obtained by writing the equations involved in linear forms of the 
type shown below is well within the accuracy of the instrument. 
* Published by permission of Humble Oil and Refining Company. 
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A pair of charts can readily be made for each instrument when used in 
72° azimuth settings and another pair for 120° azimuth settings. The 
charts appended were made for field use for a Bamberg instrument of the 
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Humble Company equipped with Blau wires and are for the latter type of 
settings. For this particular instrument the equations in this case are: 


c. 
Us 
GC. 
Usy 


II 


1.59[(As’ — As’) — (4s — A:)], 


= 2.93[(As’ — As’) + (As — As)], 


— 2 


2.75[(As’ + As’) — (As + As) ], 


mm & 


2.56[(As3' _ A.’) + (A3 + A») |. 


The first pair of equations is charted on B1 and the second on B2. An 
examination of the charts will show at once the ease and accuracy with which 


they are read. 


In accordance with the usual notation we have 


A. = "i — No, k 


1, 


ms 5 




















CHARTS FOR TORSION BALANCE READINGS 


where n; is the reading of the instrument in its kth position, and 


No = 4(m, + no + nz). 


The unprimed letters refer to the quantities bearing on the one beam, and 
the primed to those on the other beam. 

The originals of the charts are 20’’ X 20”’ with the horizontal and vertical 
scales running from —10 to +10 units at intervals of 0.1, a range well within 
that usually necessary. 


It has been found that nomographs are not quite as convenient as these 
charts for field use. 
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THE EFFECTS OF HEAT TREATMENT ON FINE 
METALLIC SUSPENSIONS 
By N. N. ZrrBew 
DEPARTMENT OF Puysics, THe Rice INstTITUTE 


(Received February 15, 1932) 


ABSTRACT 
When a suspended system is supported by a fine wire the equilibrium position 
usually changes slowly for a long time after the load is applied. The equilibrium posi- 
tion also changes with temperature. It is found that both of these disturbing factors 
can be eliminated by a suitable heat treatment of the wire. Observations have been 
made on tungsten and platinum-iridium wires of sizes suitable for use in the Eétvos 
torsion balance. 


INTRODUCTION 


ANY types of physical apparatus involve the use of a rotating system 

suspended by a thin filament of metal or other material. Numerous ir- 
regularities in the behavior of an arrangement of this sort are often observed, 
such as changes in the equilibrium position with time and temperature, 
changes in the torsion constant of the suspending filament and failure of the 
suspended system to return to its initial position after a displacement. These 
irregularities are especially serious in the Eétvos torsion balance, in which 
the suspension is very thin and is heavily loaded. Investigations of various 
sorts on the behavior of suspension wires have been published by Shaw and 
Lancaster-Jones,' Boys,’ Threlfall,’ Kénigsbergert and others. The present 
work is devoted to an attempt to remove by suitable heat treatment the 
changes in equilibrium position with time and temperature. On this particu- 
lar phase of the subject very little has been published. 


APPARATUS AND PROCEDURE 


Observations have been made on several tungsten and platinum-iridium 
wires of diameters and lengths suitable for torsion balance use. The furnace 
used to anneal the wires and for the time drift and temperature coefficient 
determinations was a brass tube wound with an electrical heating coil and 
well insulated with asbestos. A glass window inserted in the side allowed the 
angular position of the suspended system to be observed with a lamp and 
scale, and a thermometer inserted at the top allowed the temperature to be 
read. In all the work the scale distance was 120 cm and the scale divisions 
millimeters. 


1H. Shaw and E. Lancaster-Jones, Proc. Phys. Soc. London 35, 151 (1923). 
2 V.C. Boys, Phil. Mag. 23, 489 (1887). 

8 Richard Threlfall, Phil. Mag. 30, 99 (1880). 

* J. Kénigsberger, Zeits. f. Physik 40, 729 (1927). 
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The suspended system used was a cylindrical brass weight with a mirror 
attached since this design would minimize the effects of convection currents. 
The weight of the suspended system was about half that required to break 
the wire. The upper end of the suspension was fastened to a brass rod which 
extended through the top of the furnace and was clamped outside so that a 
twisting or slipping of the furnace would not affect the position of the sus- 
pended system. 

In the case of tungsten wires carbon dioxide gas was kept flowing through 
the furnace during all annealings at temperatures higher than 200°C. In 
order to anneal at temperatures higher than 500°C the wire was suspended 
in a glass tube with a weight of 6 to 10 grams at the lower end dipping ina 
copper sulfate solution so that a current could be passed through the wire. 

The temperature coefficient, which is defined as the change in the equili- 
brium position per unit temperature change, was measured by heating the 
furnace about 20°C above room temperature and allowing it to cool slowly; 
the scale reading being noted at the beginning and at various temperatures 
during the cooling. Then a graph was plotted showing the scale reading as a 
function of the temperature. Since the relation between the scale reading and 
the temperature is linear, the temperature coefficient can be easily obtained 
from the graph. It was not possible to make readings during the heating be- 
cause convection currents disturbed the suspended system in an irregular 
way and the alternating magnetic field of the heating coil exerted a torque 
on it. The possibility of convection currents affecting the position of the sus- 
pended system even during the slow cooling was investigated by making two 
sets of readings on the same wire, one set using the furnace in the usual man- 
ner, the other set by changing the temperature of the whole room very 
slowly. The corresponding readings were the same, within experimental error, 
showing that the effect of convection currents was not appreciable. For a 
typical temperature coefficient determination see Fig. 1. The coefficient is 
reckoned positive if an increasing temperature causes a clockwise rotation of 
the suspended system, as viewed from above, and negative if the rotation is 
counter-clockwise. 


EXPERIMENTAL RESULTS 


The first observations were made at room temperature on changes in the 
equilibrium position with time. In one typical case the rate of drift immedi- 
ately after loading the wire was 4.4 mm per hour. After twenty-four hours 
the motion was still easily observable but much slower, about 0.7 mm per 
hour. In some cases a change in the direction of the drift would occur. It 
seems that this time drift is caused by a slow release of strains of some sort 
in the wire and that in a sufficiently long time a stable condition would be 
reached. Just how long a time would be required is not known but Shaw and 
Lancaster-Jones! have given a case in which the drift was still going on after 
50 days. 

The next observations made at about 200°C showed that at this tempera- 
ture the untwisting process was greatly accelerated so that the stable posi- 
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tion was reached in about 7 hours instead of many days. At higher tempera- 
tures still less time was required. During this baking process the total rota- 
tion was often 180° or more. It was also found that if a wire from which the 
time drift had been removed by baking were strained beyond the elastic 
limit by a longitudinal stress, then the time drift again appeared and could 
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Fig. 1. 


be removed only by further baking. The same result followed when the brass 
cylinder attached to the lower end of the wire was twisted through about 
half a revolution and held in this twisted position for several hours. 

The time drift is thus quite easily removed. There remains, however, the 
change in equilibrium position with the temperature of the wire, as shown 
in Fig. 1. 

The next problem is then to study the effect of heat treatment on this 
temperature coefficient, which in Fig. 1 is —0.08 mm per °C. Fig. 2 shows 
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the results of such a study. The complete procedure in obtaining these re- 
sults is as follows. The wire is hung up with a load of 20 grams and baked 
at 230°C for 6 hours. This baking removes the time drift. The temperature 
coefficient is then obtained by taking the readings shown in Fig. 1 and is 
found to be —0.08 mm per °C. This value gives the first point on Fig. 2. The 
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wire is then baked at 310°C for 2.5 hours. The temperature coefficient is then 
again measured between 20°C and 50°C and found to be —0.09 mm per °C. 
The process is then continued with successively higher baking temperatures. 
Fig. 2 shows that the temperature coefficient between 20° and 50°C becomes 
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zero after baking at a certain critical temperature of about 425°C and takes 
on increasing positive values after baking at higher temperatures. Fig. 3 
shows similar and somewhat more consistent results obtained with another 
length of the same wire. The critical temperature is about the same in the 

















6 
4 
a 
2 3 
im) 
° 
ol _,___Anneoling Temp °C daa 
200 400 600 800 1000 
3 
~Z 2 
ra) 
Material Tungsten 
-4 Diameter OO2mm 
Length 50cm 
Load 2o0 qm 
~§ Scale Dist 120 cm 











Fig. 4. 


two cases. Fig. 4 shows that for tungsten the critical temperature is consider- 
ably higher. These latter results are subject to considerable error because 
temperatures higher than 500°C were estimated roughly by the color of the 
hot wire. 
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It has been found that the duration of time during which the wire is 
baked at a certain temperature is of little importance so far as the effect on 
the temperature coefficient between 20° and 50°C is concerned. In the results 
shown, this time varied from 3 minutes to 7 hours. It has also been found 
that after a wire has once been heated to a certain temperature, subsequent 
bakings at lower temperatures have no effect. That is, the curves shown can 
be traversed in one direction only, and after a wire has once been carried 
above the critical temperature no later heat treatment can diminish its tem- 
perature coefficient. Finally, measurements made on the two halves of the 
same wire show that the temperature coefficient is uniformly distributed 
along the wire and is not due to any local irregularity. 

We have no explanation for the peculiar behavior shown in the last three 
figures. However, incomplete results on a fourth sample from an entirely 
different source show the same sort of behavior, and we believe that the ef- 
fect, whatever its cause, is genuine. Suitable precautions have been taken to 
avoid the effect of convection currents and other disturbing factors. The fact 
that in all cases so far observed, involving wires from three sources, the tem- 
perature coefficient is initially in the same direction is also difficult to explain 
and may be merely a coincidence. 

In conclusion, the writer wishes to express his thanks to Professor H. A. 
Wilson, under whose direction this work was done. 
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ON THE CORRELATION OF ISOGEOTHERMAL 
SURFACES WITH THE ROCK STRATA 


By C. E. Van OrstrRANp! 


UNITED STATES DEPARTMENT OF THE INTERIOR 
(Received February 1, 1932) 


ABSTRACT 


An instance of regional variation in Oklahoma and two cases of local variation, 
one at Long Beach, California, and the other at Salt Creek, Wyoming, have been 
selected for consideration from a large number of geothermal surveys conducted dur- 
ing the past few years by the U. S. Geological Survey and the American Petroleum 
Institute. The causes of local and regional variations are unknown. Possible explana- 
tions such as radioactivity, proximity to crystalline rocks, and transfer of heat along 
the strata are given careful consideration in attempting to explain the observed rela- 
tions between the strata and the isogeothermal surfaces. 


EGARDLESS of the source of the earth’s internal heat, the mass of the 

’ earth is so large, and the time during which it has been losing heat is so 
long, that whatever the law of increase of temperature with depth may be, 
the theoretical temperature distribution to depths of one or two miles should 
be, toa very high dégree of precision, a straight line; but, in a recent summary 
of the data from 400 deep wells, distributed in 18 states, only 5 percent of the 
depth-temperature curves could be classified as linear; 36 percent were con- 
cave to the depth axis; and 59 percent were convex to the depth axis. A curve 
of the latter type is shown in Fig. 1. The causes of these anomalies are not 
definitely known. Diminution of thermal conductivity with depth resulting 
from a corresponding diminution in the moisture content of the rocks may 
be a partial explanation of the observed convexity of the curves; and the in- 
creased conductivity of the more dense crystalline or basaltic rocks beneath 
thin sediments is a possible explanation of the concave curves. The curves 
tend to be more or less of the same type in each individual field. Not with- 
standing these anomalies, however, the curves obtained under ideal condi- 
tions are exceptionally smooth and uniform, and the resulting isothermal 
surfaces determined from them are likewise smooth and uniform. 


REGIONAL VARIATION OF ISOGEOTHERMAL SURFACES 


The general trend of the 100°F isogeothermal surface? extending over a 
distance of 100 miles from Okemah to Oklahoma City, Oklahoma, is shown 
in Fig. 2. That these surfaces are closely related to the strata is evidenced by 


1 Published by permission of the Director of the U. S. Geological Survey. 

2 John A. McCutchin, Determination of geothermal gradients in Oklahoma. Bull. Am. 
Assoc. Petroleum Geologists 14, 535-555 (1930); Determination of geothermal gradients in oil 
fields located on anticlinal structures in Oklahoma. Am. Petroleum Institute, Production Bull. 


No. 205, 19-29 (1930). 
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the fact that the crystalline rocks are found at a depth of about 2,500 feet at 
Okemah and at about 10,000 feet at Oklahoma City. The observations show 
that the isogeothermal surfaces tend to parallel the granitic bed at the base 
of the sedimentary strata. The irregularities which appear at several points 
along the general trend are apparently due to local uplifts on some of which 
oil fields are located. 


LocaL VARIATIONS OF ISOGEOTHERMAL SURFACES 


The Salt Creek and Long Beach oil domes, Figs. 3 and 4, are excellent ex- 
amples of local uplift in which the isogeothermal surfaces stand in a definite 
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Fig. 3. Salt Creek Dome, Natrona County, Wyoming. 


relation to the strata. In Fig. 3, the closed broken lines represent elevations on 
the “Second Wall Creek oil sand”. Faults are represented by broken lines 
intersecting the contours. The small dots or disks represent well locations 
and the accompanying numbers represent the depths at which a temperature 
of 80°F is reached. The heavy continuous lines represent contours on the 
80°F isogeothermal surface. Fig. 3 shows clearly that the highest point on the 
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isothermal surface corresponds very closely to the highest point on the 
“Second Wall Creek sand”. Extending laterally from this point, the 100 and 
120°F isothermal surfaces tend to parallel the surface of the sand as shown in 
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Fig. 4. Long Beach Dome, Los Angeles County, California. 


Fig. 5.5 In general, the isothermal surfaces are flatter than the surfaces of the 
strata and the flatness increases with diminishing depth until a point is 
reached near the surface of the ground at which the isothermal surfaces 
tend more or less to conform to surface topography. As the effects of surface 
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Fig. 5. Cross-section of Salt Creek Dome, Wyoming. 


topography, however, are practically negligible in this field, we must look to 
other causes for an explanation of the marked rise of the thermal surfaces as 
they pass over the dome. 

The Long Beach dome’ presents a problem in which surface topography is 


*W. T. Thom, Jr., Relation of earth temperatures to burried hills and anticlinal folds. 
Economic Geology 20, 524-530 (1925). 

4 A. J. Carlson, Geothermal variations in oil fields of Los Angeles Basin, California. Bull. 
Am. Assoc. Petroleum Geologists 14, 997-1011 (1930); Geothermal conditions in oil producing 
areas of California. Am. Petroleum Institute, Production Bull. No. 205, 109-139 (1930). 
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an important factor. As shown in Fig. 4, the isothermal surfaces are excep- 
tionally smooth and uniform and follow the top of the upper oil producing 
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Fig. 6. Depth-temperature curve, deep well No. 2, Tucker and 
Johnson lease, Long Beach, California. 


strata, the Lower Brown Zone, with marvelous precision. That uniformity in 
these surfaces was to be expected is suggested by the depth-temperature 
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Fig. 7. Longitudinal section of Long Beach oil field. (See 4, first reference.) 


curve of deep well No. 2, Fig. 6. This curve like many others in the Long 
Beach field, shows a slight rise in temperature at a depth of about 1000 feet, 
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beyond which the curve is slightly convex towards the depth axis. Another 
point of interest is the fact that the excess of annual mean soil temperature 
just beneath the surface of the ground over annual mean air temperature 
just above the surface of the ground (a-b, Fig. 6), amounts to more than 
7°F. This large discontinuity in the depth temperature curves at the surface 
of the ground is typical of oil field areas in Southern California. 


CORRECTION FOR SURFACE TOPOGRAPHY 


The longitudinal section of the Long Beach field, Fig. 7, shows at once 
that the total rise in the isothermal surfaces as they pass over the dome and 
beneath the hill can not be explained on the basis of normal cooling; for, the 
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Fig. 8. Transverse section of Long Beach Dome showing computed isotherms. 


fact that the thermal surfaces are steeper and show a greater rise from the 
lowest to the highest points than the topographic surface above them proves 
conclusively that we are not dealing with a normal flow of heat towards the 
surface of the earth. In order to form an estimate of the discrepancy resulting 
from this abnormality, use will be made of the equations given by Lees.’ 

Referring to Fig. 8, let the axis of x be taken at sea level, with the origin 0 
directly beneath the apex of the symmetrical hill. Elevations above sea level 
are negative; below, positive. As the hill or mountain is assumed to be a ridge 
which is very long in comparison with its cross-section, the problem is re- 
duced at once to the flow of heat in two dimensions. The fundamental equa- 
tion to be solved is 

Oy =v w 


oe m 


’ Charles H. Lees, On the shapes of the isogeotherms under mountain ranges in radioactive 
districts. Proc. Roy. Soc. London A83, 339-346 (1910). 
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subject to the condition 


9 


k— = hv — 0 ) 
(v — Ua) (2) 


which represents the loss of heat from the surface of the earth at any point on 
the plain or mountain slope. 

In these equations, v=temperature at point x, z; v= temperature of air 
in contact with the surface of the ground at the given point; w=number of 
calories of heat generated per second per cubic centimeter by radioactive 
processes; k= coefficient of thermal conductivity of the rocks; h =coefficient 
of emissivity; » = numerical magnitude measured inward along the normal to 
the surface of the ground. 

Now let us put, also, v7) =annual mean temperature of soil just beneath 
the plain surface of the ground; a=temperature gradient beneath the plain; 
a’ = temperature gradient in the air along the mountain slope. 


(1 = (a—a’)/a 3) 


ce = (v — w)/a 


IIT=elevation of top of hill referred to level of plain; )=half-breadth of hill 
at elevation h nearly one-half the height of the hill; d=distance from apex 
of hill to point on contour at elevation h above the plain. (See Fig. 9). 


B = w/2k(a — a’) 
h (1+ 8h 
+h (4) 
H (1+ 8H) 
The quantity 6 is dependent on radioactivity. Eq. (4) is an arbitrary as- 
sumption introduced for convenience in handling the equations. 
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Fig. 9. Transverse section of imaginary hill or mountain. 


Without going through the necessary steps of integration of Eqs. (1) and 
(2) and the subsequent transformation of the integrals in order to take into 
account the cross-section of the mountain, we write at once the equation to 
the contour of the transverse section, 


ee ee ee ee 4) 0 (5) 
2(1 — Bz ~~ 5 
o ++ 8+ o 
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and the corresponding equation to the isogeotherms, 
we? » Hd(i + BH)(z + H + d) 


1— % = az—-—+(a-a 
2k x*+(:+ H+ d)? 
As Eqs. (5) and (6) are cubics in 2, it will be more convenient for purposes 
of calculation to put them in the respective forms: 





(6) 





Hd 
v--b+Hta[e+n+a+—| (7) 
c,d 
v= — [s+H+d] | +H+d)+ | (8) 
Z—- Ce 


in which the terms dependent on radioactivity have been neglected. 
Putting x =0 in (8), 


2+ (H+d—c)s — (H+ d)a+ Hd =0 (9) 


and assigning particular values to cz, which according to Eq. (3) represents the 
depth beneath the plain at which a particular temperature v is reached, Eq. 





TABLE I. Depth of isogeotherms beneath level of plain and rise of isogeotherms beneath apex of hill. 
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0 | 89.2 | 293 38. 70 | 1149. 





0 0.0 2 9 2| 3770 | 24.0| 79 
a. 2 32.8 108 | 81.8 | 269 | 41.7 75 | 1231.2 | 4039 | 22.8 | ve 
y Sie 4 | 65.7 215 | 75.8 | 249 44.4 80 | 1313.3 | 4309 | 21.8 71 
: 6 98.5 323 | 70.7 | 232 47.2} 85 | 1395.4 | 4578 | 20.9 68 
4.4 8 131.3 431 | 66.3 | 218 50.0 90 | 1477.5 4847 | 20.0 | 66 
5.6} 10 164.2 | 539 | 62.5 | 205 | 55.6 | 100 | 1641.6] 5386|18.5| 61 
6.7 12 197.0 | 646 | 59.1 194 61.1 | 110 | 1805.8 | 5925 17.1 56 
7.8 14 | 229.8 | 754 | 56.1 184 66.7 | 120 | 1970.0 6463 | 16.0 | 52 
8.9 16 | 262.7 | 862 | 53.5 | 175 | 72.2 | 130 | 2134.2 | 7002 | 15.0 49 
10.0 18 295.5 | 969 | 51.1 | 167 | 77.8 | 140 | 2298.3 7540 | 14.1 | 46 

| i | 
11.1 |; 20 | 328.3 | 1077 | 48.9 | 160 | 83.3 | 150 | 2462.5 | 8079 | 13.3 | 44 
13.9 | 25 | 410.4 | 1346 | 44.1 | 145 | 88.9 | 160 | 2626.6 | 8618 | 12.6 41 
16.7| 30 | 492.5 | 1616 | 40.3 | 132| 94.4 | 170 | 2790.8| 9156 | 12.0| 39 
19.4 | 35 | 574.6 | 1885 | 37.1 | 122 | 100.0 | 180 | 2955.0 | 9695 | 11.4 | 38 
22.2 | 40 | 656.7 | 2154 | 34.4 | 113 | 105.6 | 190 | 3119.1 | 10233 | 10.9 | 36 
| | 
25.0 | 45 | 738.7 | 2424 | 32.1 | 105 | 111.1 | 200 | 3283.3 | 10772 | 10.5 | 34 
27.8 | 50 | 820.8 | 2693 | 30.0 98 | 116.7 | 210 | 3447.5 | 11311 | 10.0 | 33 
30.6} 55 902.9 | 2962 | 28.2 93 | 122.2 | 220 | 3611.6 | 11849 9.6 | 32 
33.3 60 985.0 | 3232 | 26.7 87 | 127.8 | 230 | 3775.8 | 12388 9.3 | 30 
36.1 65 | 1067 .1 3501 | 25.2 83 | 133.3 | 240 | 3940.0 | 12926 8.9 29 
a=0.000 338 412°C per cm (=1°C in 29.55 meters=1°F in 53.86 feet). 
a’ =0.000 05°C per cm (=1°C in 200 meters=1°F in 364.5 feet).® 
H =360 feet = 109.73 meters. 


b=3.85H =1385 feet =422.15 meters. 
d= 1396.65 feet =425.70 meters. H/d=0.257 760. 
c;, =0.852 251 co=1641.654 cm per degree Fahrenheit. 


6 W. J. Humphreys, Physics of the Air. Sec. ed., p. 40 (1929). 
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(9) gives the elevation beneath the apex of the hill of that particular isotherm 
which is asymptotic to c. The difference between z and C2 is equal to the rise 
of the isothermal surface above its asymptotic plane. Thus in Fig. 8, oa rep- 
resents the rise of that isothermal surface which is just beneath the surface 
of the plain, (cg=0), and which has ox extended for an asymptote; bb’ repre- 
sents the rise of that surface whose asymptote is perpendicular to the axis 
of the hill at b’; and so on. These values are tabulated in Table I. The table 
applies to any other hill of the same cross-section and the same geothermal 
gradients (a and a’) beneath the plain and along the mountain slope. To 
determine the temperature represented by any isotherm, it is merely neces- 
sary to add the tabular value in the table, or Fig. 8, to the annual mean 
temperature of the ground just beneath the surface of the plain. 

In order to compare the theoretical rise of the isogeothermal surfaces 
with the observed, the gradient beneath the plain was assumed to be, a@ 
=(0.000,338,412, corresponding to a rate of temperature increase of 1°C in 
29.55 meters, or 1°F in 53.86 feet. This value is an approximate average of 
the gradients for 42 wells in the Long Beach field. From the same wells, the 
average annual mean temperature just beneath the surface of the ground was 
found to be, 7=71.6+0.2°F. Hence, the 80°F isotherm corresponds to a 
temperature difference in Table I of 8.4°F. Interpolating, we find the rise of 
the isotherm to be 215 feet. The last two columns of Table II were obtained 
by a repetition of this process, using for the last column a table similar to 
Table I. In column 4, the level of the plain was assumed to coincide with sea 
level; in the last column of the table, it is taken to be 40 feet above sea level. 
The data for the transverse section of the hill, Fig. 8, were kindly forwarded 
to me by Dr. A. J. Carlson. The observed differences in elevation of the sur- 
face of the ground and the isotherms, column 3, Table II, are taken from 
his paper on “Geothermal variation in the oil fields of Los Angeles Basin, 
California.” 

Comparison of column 3 of Table II with the last two columns shows 
that the observed heights are somewhat more than twice the computed 


TABLE II. Rise of isogeotherms. 























Differences in elevation 
Depth 
Isogeotherms beneath _ Assumed | Assumed 
F = Observed Computed | Computed 
feet feet feet 
296 | 360 | 320 
80 453 492 215 | 189 
90 991 348 166 144 
100 1530 274 | 136 | 118 
110 2068 223 116 100 
120 2607 | 182 | 100 | 86 
130 3146 | 200 | 89 | 76 
140 3684 141 80 68 
150 4223 | be | 72 | 62 
Lower Brown | 
Zone — | 1100 — — 
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heights. In making this comparison, I have attempted not to underestimate 
the computed heights. For example, the assumed elevations of the top of the 
hill above the level of the plain, namely 360 and 320 feet, instead of 296 feet 
given by Dr. Carlson, gives a rise of isothermal surface which is too high. 
Likewise, the substitution of a long ridge for a hill tends to make the com- 
puted values too high. An important possible exception, however, should be 
noted. Let us assume that the gradient beneath the plain is the steepest in 
the field corresponding to a rate of 1°F in 48.5 feet instead 1°F in 53.86 feet 
as used in constructing Table I. Hence we have 


a = 0.000 375 811 c, = 0.866 954 


and substituting in Eq. (9), putting c.=0 in order to represent the 0°F iso- 
therm, we find 299 feet instead of 293 feet as recorded in Table I. This dif- 
ference of 6 feet is rather insignificant in comparison with the large difference 
of more than 100 feet for which we require an explanation. Using the value 
a’ =0.00006 instead of a’ =0.00005, the corresponding rise is 288 feet instead 
of 299 feet. This shows that the value of a’ is of importance in determining 
the computed heights. Taken as a whole, however, the evidence substantiates 
the conclusion that the differences in elevation between the highest and low- 
est points on the observed isotherms are abnormally large. 


COMMENTS ON THE DATA OF OBSERVATION 


The tests in the Salt Creek field were made under conditions such that a 
close approximation to true rock temperatures is to be expected. Values of 
a in the equation 


y=at+bdx 


in which y=temperature at depth x; a=annual mean temperature just be- 
neath the surface of the ground; )=gradient, vary from 46.9 to 52.3°F. The 
mean from 19 wells is 49.8+0.2°F. The probable error r of an observation of 
weight unity is 1.0°F. Values of the excess of soil temperature over air tem- 
perature (a—b, Figs. 1 and 6) range from 0.2 to 5.6°F. The mean is 
3.10+0.2°F. r= +1.0°F. 

At Long Beach, the temperatures just beneath the surface of the ground 
from 42 wells fall between the limits 68.4 and 77.2 with a mean value of 
71.6+0.2°F. r= +1.4°F. The maximum value of the excess of soil tempera- 
ture over air temperature is 14.6, the minimum, 5.8°F. The mean is 9.0 
+0.2°F. r= +1.4°F. 

The large range in the preceding values suggests the possibility that the 
temperatures in some of the wells that had been pumped are too high. The 
resulting error is a maximum at or near the surface of the ground and gradu- 
ally diminishes with more or less regularity to zero at the source of the rising 
column of oil or water. An obvious effect of this abnormality is to elevate the 
isotherms in amount roughly proportional to their elevation above the fluid 
bearing sand, and diminish the gradients or increase the reciprocal gradients 
in the same relative proportion. To ascertain the possible magnitude of this 
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TABLE III. Mean values of a, b, and 1/b. 












































100-1000 feet | 100-2000 feet 
No. | | | | | | No. | | 
of | a r | ro | b 7 X108r9X105, 1/b | of b ir X 10579105) 1/b 
Obs. | | | Obs. | 
5 | 72.8} 1.2 | 0.5 0.01781 | 58 26 | 56.1! 5 0.01776 | 44 20 | 56.3 
35 | 71.3 | 1.4 | 0.2 | .01885 9 2 | 53.0 | 35 | .01878 | 59 10 | 53.2 
; | | | 
10 | 73.4 | 1.6 | 0.5 |0.01804 | 87 28 | 55.4 | 10 0.01804 | 63 20 | 55.4 
30 | 70.91 1.0} 0.2 | .01895 | 88 16 | 52.8 | 30 | 01885 | 55 10 | 53.0 
| | 

10 | 73.4 | 1.6 | 0.5 |0.01804 | 87 28 | 55.4 10 \0.01804 | 63 20 | 55.4 
10 | 70.9 | 1.0 | 0.3 | .01919 81 26 | 52.1| 9| .01923 | 34 11 | 52.0 
10 | 71.2 | 1.1 | 0.3 | .01855 | 108 34 | 53.9 | 10 | .01849 | 69 22 | 54.1 
10 | 70.4 | 1.1 | 0.3 | .01912 | 74 23 | 52.3 | 9 | .01882 | 38 13 | 53.1 

















| 





source of error, Table III has been computed. The table is based on the data 
from 40 wells, two of the poorest values being omitted. In the first two lines, 
the mean values for 5 wells highest on the structure are compared with the 
means from the remaining 35 wells. Similarly, in the third and fourth lines, 
the means for 10 wells highest on the structure are compared with the means 
from the remaining 30 wells, while in the final tabulation the wells have been 
grouped in four successive groups as regards elevation on structure, the first 


TABLE IV. Gradients and reciprocal gradients (H = 360 feet; and H =320 feet). 





















































In air on slope Beneath adjacent Beneath summit 
of hill plain of hill 
Gradients Reciprocal Gradients Reciprocal Gradients Reciprocal 
(a’) gradients (a) gradients (as) gradients 
°C per Feet | Meters | °C per Feet | Meters °C per Feet | Meters 
cm per °F | per °C cm per °F | per °C cm per °F | per a 
0.00006 303 .8 166.7 0.00036 50.6 27.8 0.00028 27 64.5 35.4 
0.00006 | 303.8} 166.7 | 0.00030 | 60.8 | 33.3 | 0.00023 81 76.5 | 42.0 
0.00005 | 364.5 | 200.0 | 0.00036 | 50.6| 27.8 | 0.00028 01} 65.1 | 35.7 
0.00005 | 364.5 | 200.0 | 0.00030 | 60.8 | 33.3 | 0.00023 56| 77.4} 42.5 
0.00006 303 .8 166.7 0.00036 50.6 27.8 | 0.00028 64 63 .6 34.9 
0.00006 303 .8 166.7 0.00030 60.8 33.8 0.00024 11 75.6 41.5 
0.00005 | 364.5 | 200.0 | 0.00036 | 50.6 | 27.8 | 0.00028 40 | 64.2 | 35.2 
0.00005 | 364.5 | 200.0 | 0.00030 | 60.8 | 33.3 | 0.00023 87 76.4 | 41.9 























group as heretofore being the highest. The table shows clearly that the values 
of a and 1/0 in the first group of each set are always the highest. However, 
minor irregularities of this kind can not always be attributed to pumping 

operations. In general, the depth-temperature curves become more irregular 
: as the surface of the ground is approached. It is for this reason that a map 
showing the isothermal surfaces is much more reliable than a map showing 
the distribution of reciprocal gradients. In addition to the preceding irregu- 
larities, the variation of the annual mean temperature of the air over the 
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surface of the hill must be taken into account. Thus, the value of a for six 
wells on the south side of the hill is 73.0+0.8; the same for 15 wells on the 
north side is 71.4+0.3, making a possible temperature difference on the two 
slopes of 1.6+0.9°F. Other observers have likewise found maximum tem- 
peratures on southern slopes.’ 

As a further explanation of the slight increment in the observed values 
of (1/8) as the crest of the hill is approached, a comparison must be made 
between the assumed gradient (a) beneath the plain and the theoretical 
gradient (a,) beneath the summit of the hill. Using the value of H/d given 
in Table I and Fig. 8, and making the appropriate substitution in the equa- 
tion 


a,=a— H(a—a’)/d 


we obtain the values of a, tabulated in the first part of Table IV. In the sec- 
ond part of the table, (@,) is given for H=320 feet, H/d=0.245 240. 

Interpolating for our assumed reciprocal gradient beneath the plain, 
1°F in 53.86 feet, we find the value 1°F in 68.3 feet at the summit of the hill. 
That is, the computed variation from the plain to the summit of the hill is 
four or five times the differences, 3.1 and 2.3, shown between the highest and 
lowest groups on the structure in the last set of values in Table III. In one 
way only does it seem possible to bring theory and observation into agree- 
ment. Substitution in the equation, 





Ov Hd 
Oz (c+ H+ d)? 


which represents the gradient at any elevation (z) beneath the apex of the 
hill, we find for 


a = 0.00033 8412(1° F in 53.86 feet), a’ = 0.00005, H = 360 feet, 


the values, 


1/6 =1°F in 62.5 feet z=0 meters 
= ” 55.9 ” 2=500 od 
= ” 54.8 ” z= 1000 


At z=o0, 360 feet beneath the top of the hill, the theoretical difference of 
8.6 feet between the assumed and computed values of (1/b) greatly exceeds 
the observed differences, 3.1 and 2.3 feet, Table III; but at z=500 meters, 
corresponding to a depth of 2000 feet beneath the top of the hill, the observed 
differences slightly exceed the computed value of 2.0 feet. At greater depths, 
Eq. (10) shows that the gradients beneath the plain and the hill approach 
coincidence. 

To sum up the evidence at Long Beach, the elevation of the isotherms 
with reference to the topographic surface, and the small variation of gradient 
over the structure in comparison with the theoretical variation proves con- 


7 Edith M. Fitton, Soil temperatures in the United States. Monthly Weather Review 59, 
9 (1931). 
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clusively that the temperature distribution in this field is highly abnormal. 
The only exception appears to be the possible agreement of the observed and 
computed gradients at depths of about 2000 feet. 

At Salt Creek, the evidence of variation of temperature with structure 
is conclusive. At Long Beach, the evidence is confused on account of the fact 
that in some of the wells the temperatures at the higher levels were increased 
by pumping operations while in a few of the remaining wells, the tempera- 
tures may have been lowered as a result of the wells being on a vacuum. 
However, proceeding on the assumption that we have a normal distribution of 
temperatures, and taking into account the fact that the abnormal tempera- 
tures in some of the wells have tended to increase the values of (1/0), it 
follows that the observed variation in the values of (1/b) as we pass from the 
plain to the summit of the hill, should have exceeded slightly the theoretical 
variation; but, just the converse is true—the observed variation is only a 
small fraction of the theoretical variation. Hence we are driven to the con- 
clusion that abnormal temperature conditions in this field have tended to 
neutralize the effects of surface topography. 

















Fig. 10. Long Beach Dome. The numbers opposite well locations 
represent feet per degree Fahrenheit, 100-2000 feet. 


CAUSES OF TEMPERATURE VARIATIONS 


Depth to the crystalline rocks and flow of heat along the strata appear to 
be possible explanations of the close relationship that apparently exists be- 
tween the isotherms and the general trend of the strata in central Oklahoma 
(See Fig. 2). 

A similar explanation may apply at Salt Creek, which is an exceptionally 
well defined local uplift. Steep dips occur on the western side of the field 
(See Figs. 3 and 5) and very hard granite was found by The Midwest Re- 
fining Company at a depth of 5420 feet near the top of the dome in the S.W. 
1/4 N.W. 1/4, sec. 35, T. 40 N., R. 79 W. This evidence suggests the pos- 
sibility that the depth to the granite is least on top of the dome, and that the 
direct contact of the granite with the hot crystalline rocks at great depths 
serves to maintain lines of flow from an almost inexhaustible heat source to 
the top of the dome. Additional evidence in substantiation of this conclusion 
is to be found in the fact that the values of 1/b vary from 1°F in about 29 
feet on top of the dome to 1°F in about 56 feet in Sec. 15, T. 40 N., R. 79 
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W. on the outer edge of the dome. Not only is the great range in values 
of 1/b in agreement with our hypothesis, but the further fact that abnor- 
mally high temperatures are found immediately above the granitic mass 
leaves little to be desired in explanation of the observed distribution of tem- 
peratures. 

On account of the great depth to the crystalline rocks at Long Beach, 
possibly 20,000 feet, it seems unlikely that depth to the crystalline rocks, 
which appears to be an important factor at Salt Creek and in Central Okla- 
homa, will provide an explanation of the temperature variations found in 
this field. Steeply dipping strata are present, but about the only other pos- 
sibilities to which we can appeal are recent erosion, abnormal thermal con- 
ductivities and generation of heat in the strata. I shall not here attempt to 
discuss the merits of these possibilities. 

To estimate the effects of radioactivity on the isogeotherms at Long 
Beach, let us expand Eq. (4) into the form 


H 8 
w= S|i+oat- | 
2 wy) 


- 


and then putting 


| 


H = 360 feet = 10973 cms w= 1.5 X 10° 


0.00036 a’ = 0.00006 k = 0.004 


a 


we have 
HH | w 
h = — = 5486 cms 6 = ————— = 6.25 X 1078 
2 2k(a — a’) 
and the new value of / is found to be 5488 cms an increase of only 2 cms. : 
Using Joly’s* recent value, w=42.110-“, we have 


B=1.75X* 10-7? hk = 5491 cms 


an increase of 5 cms. The resulting change in d (Fig. 9) is obviously negligi- 
ble, consequently Eq. (6), which contains (d), remains unchanged. 

Another method of procedure consists in computing the rise in the 0°F 
isotherm beneath the apex of the hill. Putting «=o and z=o in (6) and using 
the same constants as before, including Joly’s value of w, we find 





+A ‘(1 + BH : 
i v= (a — a’)( BH)- Hi 
1+— 

d 


2.617 + 0.005 = 2.622°C 


from which it again appears that the term due to radioactivity is negligible 
in comparison with the observed variations in temperature across the struc- 
ture. 


8 John Joly. The surface-history of the earth. Second edition, p. 60 (1930). The Clarendon 
Press, Oxford. 
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SoME GENERAL SUGGESTIONS 


Our observations are not yet sufficiently accurate nor sufficient in num- 
ber to permit making definite generalizations. The chief source of inaccuracy 
is unstable temperature equilibrium in wells. Notwithstanding these irregu- 
larities, however, the observations show conclusively that the isogeothermal 
surfaces rise in passing over a large number of salt domes and anticlinal struc- 
tures. There is no well established evidence of depressed isogeotherms in any 
oil field, but in a considerable number of fields, the rise in temperature is not 
very marked. In general, the observations suggest the conclusion that a rela- 
tion exists between the high temperatures and positive gravity anomalies 
that prevail over the crests of many of the oil field anticlines. 

The causes of temperature variations given in the preceding paragraphs 
are to be regarded merely as preliminary suggestions. The complete explana- 
tion will undoubtedly involve isostasy, geological history, and the sources 
of the earth’s internal heat. 
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ABSTRACT 


The Michigan College of Mining and Technology, in cooperation with the 
Calumet and Hecla Copper Company and the author, is carrying out a program of 
temperature measurements in the deep copper mines of Northern Michigan, extending 
the previous work of Agassiz and others. Temperatures are measured with mercury 
thermometers mounted in bakelite tubes, placed in drill holes in mine workings where 
the rock has been freshly exposed, special attention being given the effects of drilling, 
blasting, and other heat conduction considerations. Present results give as the average 
gradient from the surface to 5679 feet below (temp. 95.3°F), 1°F in 108.5 feet (0.0168° 
C/meter). The gradient is more nearly uniform than has sometimes been supposed. A 
preliminary attempt has been made at calculating the previous “thermal history” of 
this region. Diffusivity of specimens of the rock measures 0.0075 c.g.s., and on this 
basis calculations of theoretical temperature-depth curves have been made for 25 
different assumptions of previous temperature conditions, and compared with the 
actual curve. Results as yet are inconclusive but indicate that at least 30,000 years 
have elapsed since the last glacial epoch, a longer period than usually assumed. 


VER since the days of Kelvin’s theoretical discussion of the age of the 

earth there has been interest in measurements of the earth’s temperature 
gradient, and while radioactive discoveries have considerably diminished the 
importance of calculations on this basis there are nevertheless many other 
interests—some of them new—attaching to geothermal measurements. As 
a result a mass of thermal data has been and is being collected from deep 
mine, well and drill hole measurements in all parts of the world. The copper 
mines of the Keweenaw peninsula in northern Michigan offer a particularly 
inviting opportunity for such study because of their great depth and freedom 
from local heat sources. The pioneer measurements of Agassiz,! which re- 
sulted in an apparently unusually low gradient, are too well known to require 
comment, and Lane,? Van Orstrand’ and others have also made valuable con- 
tributions along this line. 

Certain favorable opportunities, depending partly on the increasing depth 
of the mines, have for some time called for a reopening of this study and ac- 
cordingly Professor James Fisher of the Michigan College of Mining and 
Technology and Mr. Harry Vivian, Chief Engineer of the Calumet and Hecla 





1 Alexander Agassiz, Am. Jour. Sci. (3) 50, 503 (1895): Brit. Assoc. Adv. Sci. Report of 
71st meeting, p. 65 (1901). 

2 A.C. Lane, Bull. Geol. Soc. of Am. 34, 703 (1923). 

$C. E. Van Orstrand, Am. Jour. Sci. 15, 495 (1928). See also N. H. Darton, U. S. Geol. 
Survey, Bulletin 701, p. 50. 
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Consolidated Copper Company, are cooperating with the writer in a program 
of work which may extend over a period of years. The Michigan College is 
generously supplying the financial aid and the Calumet and Hecla Company 
the skilled technical assistance necessary to the success of the work. The 
measurements in all cases are to be guided by heat conduction considerations 
to the end that they may represent as closely as possible the actual virgin 
temperature at the spot, unaffected by mining operations. The work must 
accordingly proceed slowly, waiting frequently for favorable opportunities 
for temperature measurement; but while it has been going only about a year 
some preliminary measurements of value have been obtained as well as other 
results which may be of interest to the physicist. 


Metruop or TEMPERATURE MEASUREMENT 
General procedure 


The measurements so far have been confined to the Calumet and Hecla 
mine in which most of the Agassiz measurements were also made. A some- 
what different procedure has been adopted, however. In the Agassiz work 
thermometers were sealed for weeks at the bottom of drill holes 10 feet deep, 
in shafts or passages which had, themselves, been exposed to ventilation in 
some cases for many months. A little calculation on the basis of heat con- 
duction theory will serve to show that the results under such conditions are 
very likely to be influenced by the ventilation, and only in the most favorable 
cases can be depended on to give virgin rock temperatures. 

In the present work the measurements have all been carried out in new 
workings, i.e., drifts or other cuttings which have been advancing steadily 
a number of feet a week and which are well removed from other parts of the 
mine. At such a “temperature station” a hole is drilled a few feet back from 
the breast and in rock whose face has been exposed only a few days. The hole 
is usually 7 feet deep, but in some cases special 14 foot holes are run. Tem- 
peratures are taken with two or more thermometers located at the bottom of 
the hole. It may be remarked that, save for the water used in drilling, which 
quickly drains out, the hole has been found practically dry in all cases. 


Thermometers 


After careful consideration of the advantages and disadvantages of elec- 
trical temperature measuring instruments, mercury thermometers have been 
finally chosen for this work, at least for the type of measurements being made 
at present. They have been specially made by Henry J. Green, reading 
0—40°C in 1/10°, and frequent zero tests and comparisons with one of the 
number calibrated by the Bureau of Standards give maximum errors of the 
order of 0.04°C. They are mounted in bakelite tubes 1 inch in diameter and 14 
inches long with the thermometer bulb specially insulated thermally. This is a 
vital matter to which a great deal of attention has been given. It is necessary 
that after pulling from a drill hole the mercury thread show no change for 
something like a minute, so as to give the operator plenty of time to make his 
readings. At the same time there must be combined with this lag or “delay 
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action” a reasonably rapid action of the thermometer, once the bulb starts 
to change in temperature, as otherwise it would occupy too much time in 
coming to the temperature of the hole. 

After some calculation and much experimentation a very satisfactory 
solution of the problem has been reached, involving insulating the bulb in 
vulcanite, wax and dynamo paper with waterproof lacquer as shown in Fig. 
1. If such a thermometer is subjected to even such a radical temperature 
change as is involved in plunging into ice water, it is a full minute before the 
mercury thread shows a change of 1/10°. Having started to move, however, 
it falls several degrees a minute and comes down (asymptotically) to the 


Rokelire Vulcanite 


tube Paper 





Fig. 1. Section of mounted thermometer (Overall length, 35 cm), 


temperature of the bath—within 0.02°, say—in about half an hour. The 
final design of these thermometers has proved very satisfactory indeed for 
the work in hand. Difficulties experienced with the early ones over the se- 
paration of the mercury thread have been entireiy overcome in the later ones 
by having them made with about one atmosphere of inert gas pressure over 
the mercury. This has the further advantage of relieving the bulb of practic- 
ally all pressure, inside or outside, and doubtless tends to render the zero 
still more stable. 


Temperature readings 


As the holes are drilled in the usual manner with air-operated machines, 
with a flow of water through the drill steel, there is almost no heating effect 
due to drilling. Temperature readings taken immediately after the hole is 
completed are sometimes a degree high, but after a few hours all evidence of 
drilling heat has usually disappeared. To be on the safe side, however, little 
reliance is placed on readings made less than 24 hours after drilling, and the 
whole question of possible small errors due to heat of drilling is one on which 
considerably more study is to be expended. 

Thermometers are used in tandem groups of two or three, run into the 
hole with a stout wire, with the last one carrying a rubber gasket. The mouth 
of the hole is stopped with a wooden plug. While the thermometers reach 
equilibrium in half an hour they are left for at least an hour and then quickly 
withdrawn and read. They are then replaced in the hole, to be read again 
about every two hours until several sets have been taken. This procedure is 
repeated a couple of days later, with occasional readings continuing for sev- 
eral weeks or months. In addition, air, rock face and psychrometer readings 
are occasionally taken in the drift to furnish general information as to the 
thermal conditions at the station. 

Plots of these temperatures with time usually give an approximately 






























GEOTHERMAL GRADIENT DETERMINATIONS 157 





straight line, sloping slowly downwards, showing the gradual cooling—of the 
order of a degree or two a month—which takes place as heat is conducted 
from the rock into the slightly cooler drift. A slight extrapolation of the curve 
gives the initial (virgin) rock temperature with an accuracy believed to be 
of the order of 1/10°C. To make sure that drilling and blasting effects are 
really negligible, as any calculations which can be carried out would indicate, 
an occasional 14 foot hole has been run in at the bottom of a short cross-cut 
such as occasionally occurs in drifting. In such a hole the tests made so far in- 
dicate a steady temperature for weeks and months and go far towards remov- 
ing all doubt in one’s mind as to any transient effects. The agreement which 
has been found so far between measurements in such deep holes and the 
standard 7 foot ones is good, but much more work remains to be done along 
this line before complete confidence can be placed in the results, however 
satisfactory they appear to be as a whole. 


RESULTS 


The results to date are shown in Fig. 2. The points determined in the 
present work are seen to lie very well on a smooth curve which is almost a 
straight line. When this is continued to the surface level, making such use 
as seems justifiable of the Agassiz measurements, and especially of the point* 
indicating average surface temperature, we have the geothermal curve for 
this particular spot—insofar as it can be considered to be established by 
these measurements; it is unfortunate that there have been no recent op- 
portunities for reliable measurements in the upper mine levels. The absence 
of local heat sources in this region, as well as general geological considera- 
tions, would lead one to believe that this curve will probably be found to be 
typical of this region, but measurements in other mines will be made as op- 
portunities offer. A comparison with all the Van Orstrand drill-hole data for 
this locality must also be made. It may be noted in passing that the average 
gradient found here of 1°C in 59.5 meters is less than half as steep as that 
taken by Kelvin as the average for the whole earth, viz., 1°C in 27.76 meters. 


Theoretical interpretation of results 


A preliminary attempt has been made at the interpretation of this curve 
in terms of the previous thermal history of this region, involving particularly 
the time since the last ice age. In order to do this it is first necessary to de- 
termine the thermal diffusivity (thermal conductivity divided by the prod- 
uct of specific heat and density) of the rock material. B. O. Peirce’ made a 
number of careful tests of the conductivity of specimens of rock from the 
Calumet and Hecla mine many years ago, but unfortunately several assump- 
tions are involved when it comes to getting diffusivities from these results. 
Accordingly, a few tests on specimens brought from two regions in the mine 


‘ This is the point taken by Lane (reference 2, p. 703-705) as the average surface tempera- 
ture at Calumet. It was arrived at in several ways, all of which give nearly the same result and 
is believed to be a very good average. 

5 B. O. Peirce, Proc. Am. Acad. of Arts and Sci. 38, 652 (1903). 








158 L. R. INGERSOLL 


have been made by a method developed by the writer® some years ago. The 
specimen of rock, which must be in the form of a slab say 6 cm thick and 
25 or more cm square, is mounted in a paraffined wood frame which ther- 
mally insulates the edges, leaving only the faces exposed. A fine thermocouple 
is located in a small hole drilled into the center of the slab, everyth ng being 
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Fig. 2. Geothermal measurements in the Calumet and Hecla mines. Average temperature 
gradient, 0-5500 ft. depth, 1°F in 108.5 ft. or 0.00922°F /ft., or 1°C in 59.5 meters, or 0.0164° 
C/m. Gradient, 3500-5500 ft. depth, 1°F in 103.1 ft. or 0.00970°F /ft. 


carefully waterproofed. When such a slab has been left in a constant temper- 
ature region for 24 hours and is then suddenly plunged into stirred ice water 
so that the faces are quickly lowered to and kept at zero, the diffusivity can 
be readily calculated from the time it takes for the center temperature to 
fall half way from its original value to the zero point. Two specimens of trap 
rock tested in this way gave values for the diffusivity of 0.0077 and 0.0074 
respectively. Since the former specimen was from the vein and the latter 

® For a brief account of this method see L. R. Ingersoll and O. A. Koepp, Phys. Rev. 24, 


92 (1924). For the underlying theory see Ingersoll and Zobel, “Mathematical Theory of Heat 
Conduction” (Ginn) p. 105 ff. 
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from the surrounding trap rock, which composes the great bulk of the mater- 
ial, the weighted average is taken as 0.0075 c.g.s. units. 

With the above value, calculations’ have been made of theoretical tem- 
perature-depth curves for some 24 different assumed “thermal histories” of 
this region. The simplest assumption is that the surface temperature has been 
about as at present for 10,000 years and that previous to this the surface 
was covered with ice. In other assumptions the time has been changed and 
in still others several glacial overflows are taken account of. The periods of 
time have been varied between wide limits. The equation for calculating the 
points of these curves is 


Es ' : 
6 = =f F(a)em* 70-9). (t — Q)-8/2dd 
2h(m)? : 0 


where x is the distance below the surface, at which the temperature @ is to 
be calculated, h? the diffusivity, and F(A) the surface temperature function. 

These curves having been calculated and corrected for the general aver- 
age gradient slope, a careful comparison has been made between them and 
the actual curve as given in Fig. 2. The results are a little disappointing at 
first as the curve is so nearly a straight line anyway and so much depends 
on the justification for certain assumptions which cannot as yet be proved 
satisfactorily, e.g., that the diffusivity of the rock is unchanged for the first 
mile or more in depth. However, certain assumptions fit much better than 
others and point on the whole to a somewhat longer period—say 30,000 
years—since the last glacial overflow for this region, than is usually assumed. 

In addition to the acknowledgements at the beginning of this article, 
special thanks are due to Mr. James McNaughton, President of the Calumet 
and Hecla Company, for placing at the service of this investigation the facili- 
ties of the mines, and particularly to Messrs. H. E. Jefferson, H. S. Donald 
and R. F. Wilson for their care in making the measurements, as well as for 
many valuable suggestions. 


7 I am glad to acknowledge the assistance of the following of my students in making these 
long and arduous calculations: Miss M. C. Wolf and Messrs. R. O. Anderson, R. E. Erickson, 
R. G. Herb, W. L. Hole, R. W. Prucha and M. T. Rodine. 
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ABSTRACT 


The velocity of elastic waves in granite was determined at Quincy and Rockport, 
Massachusetts, and Westerly, Rhode Island. The waves measured were generated by 
dynamite explosions. They were recorded by portable seismographs at distances rang- 
ing from fifty feet to four thousand six hundred feet. The observed velocities for longi- 
tudinal waves were: 

Quincy......16,260+ 70 ft./sec. or 4.96+0.02 km/sec.* 

Westerly. ...16,400+120 ft./sec. or 5.00+0.04 km/sec. 

Rockport....16,670+ 40 ft./sec. or 5.08+0.01 km/sec. 

Average. ....16,530+ 90 ft./sec. or 5.04+0.03 km/sec. 

A three-component seismograph, used only at Quincy, recorded transverse waves, 
the velocity of which was 8150+90 ft/sec., or 2.48+0.03 km/sec. From the two 
velocities determined at Quincy and the density of specimens taken from the shooting 
location, 2.65 grams/cm*, values for the bulk modulus, k, compressibility, 8, rigidity,u, 
Poisson’s Ratio, s and Young’s Modulus E, were obtained as follows: k =44+1 X10" 
dynes/cm?; 8=2.28+0.05X10-" cm?/dynes; »=16.3+0.4X10'® dynes/cm?; 
o =0.333+0.005; E=43 +110" dynes/cm?. The form of the time-distance curves, 
straight lines through the origin, indicated that the waves did not penetrate deeply. 
Accordingly, the values obtained are for pressures of only a few atmospheres. The 
bearings of these results upon earlier investigations of the elastic constants of granite 
are discussed. Although direct comparisons between laboratory and field results 
are not conclusive, they indicate that the Adams and Williamson curve is incorrect 
for pressures below 2,000 megabars, and that there is no marked difference between 
dynamically and statically determined compressibilities of granite. 


I. LocaTIOoNn 


ETERMINATIONS of the velocity of elastic waves in granite were 

made in three well-known granite quarrying areas, Quincy and Rock- 
port, Massachusetts, and Westerly, Rhode Island. Within each area a zone 
of practically no topographic irregularity was selected. The continuity of the 
granite in each zone was proven by numerous quarries and outcrops. Maps 
of these zones appear in Fig. 1. The observing stations and shot locations 
shown there were mapped by plane table on a scale of 300 feet to the inch. 


Il. Mretruop 


The elastic waves were produced by the explosion of charges of dynamite 
which were always placed in contact with the granite. Sixty percent quarry 
gelatine was used. Their arrival at a point on the granite at any desired dis- 
tance from the explosion was recorded by a seismograph placed at that point. 


* The + values given in this paper are probable errors. 
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The distance travelled by the waves and the time required constituted the 
observed data. Distances were scaled from the maps reproduced in Fig. 1. 
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The time plotted against the distance yielded a time-distance graph which 
was a straight line through the origin. The velocity, the reciprocal slope of 
this line, was computed by a least square solution. 
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Fig. 1. 
Ill. INstRUMENTS 


The seismograph used at Quincy consisted of three components, two 
horizontal and one vertical, mounted on a single base. The instrument was 
always oriented so that one of the horizontal components recorded vibrations 
on a line from shot to recorder, while the other, at right angles to this, 
recorded vibrations transverse to that line. These are designated respectively 
as the longitudinal and transverse components. The longitudinal component 
had a free period of 0.77 sec., the transverse, of 0.62 sec., and the vertical, 
of 0.55 sec. Each was critically damped by a magnetic device of conventional 
design. The static magnification, effected by mechanical and optical means, 
was about one thousand. 

At Westerly and Rockport, a vertical component seismograph was used. 
It had a free period of approximately two tenths of a second and was criti- 
cally damped by a piston working in a cup of oil. The static magnification, 
also mechanical and optical, was about twelve thousand. 

Both of the instruments recorded photographically on a strip of bromide 
paper one inch wide, moving at a speed of about twenty centimeters per 
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second. A vacuum-tube driven tuning fork controlled the interruption of a 
beam of light, thereby producing timing lines, which extended the full width 
of the record, at intervals of 0.01 sec. 

The electrical circuit which was used to detonate the dynamite from the 
recording location also served for the registration on the record of the in- 
stant of the explosion. A modified polarized relay in this circuit was actuated 
both when the firing circuit was closed and when it was broken by the ex- 
plosion. A mirror attached to the moving vane of the relay and in the path 
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Fig. 2. Specimen Quincy records T, transverse component; V, vertical component; 
L, longitudinal component. 


of the beam from light-source to seismograph effected the actual registration. 
In Fig. 2, the upward jog in all seismograph lines marks the deflection caused 
by closing the firing circuit, and the downward jog a few hundredths of a 
second later, the instant of the explosion. 


IV. Frrine Positions 


At Quincy and Rockport, several abandoned water-filled quarries were 
used as firing locations. The charges were placed against the quarry walls and 
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in some cases as much as thirty feet beneath the surface of the water, and 


one 
the 


hundred feet beneath the surface of the ground. This arrangement had 
important advantages that the explosions were confined by the water and 


that the energy entered sound, unweathered granite. The efficient confine- 
ment of the explosions by the water made it possible to use relatively small 
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Fig. 3. 


rges. About five hundred pounds of dynamite sufficed for the entire in- 


vestigation. 


At Westerly, shooting positions were chosen where the granite was 


covered with several feet of glacial drift. There the charges were placed in 


holes dug down to the granite. 





164 L. D. LEET AND W. M. EWING 


V. Quincy REsULTs 

Specimen records from Quincy are reproduced in Fig. 2. The traces 
marked “7”, “V”, and “L” are the transverse, vertical, and longitudinal 
components respectively. Each record is identified by figures indicating the 
distance, shot position, recording position, and weight assigned it in the 
least square solution. For example, the first record in Fig. 2 is for a distance 
of 1050 feet. The figures 58-64—4 indicate that the shot was fired at position 
58 (see map, Fig. 1), recorded at position 64, and accorded a weight of 4 in 
the least square solution. Weights were assigned on the basis of amplitude 
of trace and general photographic quality, which determined the accuracy 
with which the time could be read. 

The interval between the instant of the explosion and the “0” timing line, 
the travel-time of the longitudinal wave on the vertical component, and the 
travel-time of the transverse wave on the transverse component are shown 
on each record. 

The time-distance graph for Quincy is shown in Fig. 3. The fact that it is 
a straight line through the origin indicates that the waves did not penetrate 
deeply. It is probable that the penetration was not more than twenty meters, 
and certain that it did not exceed 200 meters.* The lines drawn through the 
points were determined by least square solutions. Their reciprocal slopes 
give: 

Vi. =16,260+70 ft./sec. or 4.96+0.02 km/sec. 
Vr= 8,150+90 ft./sec. or 2.48+0.03 km/sec. 


where V, and V7 are the velocities of longitudinal and transverse waves 
respectively. . 

The density of specimens taken from the shooting location was 2.65 
+ 0.02 grams/cm.* 

The elastic constants of the granite can be obtained by substituting these 
values in the following equations.® 


. A+ 2y\'? 
vs = (22) . 





p 
Vr = (u/p)*” (2) 
k =X + 2/3u (3) 
= 1/k (4) 
: u(3rA + 2p) . 
«0 Se (5) 

A+u 
IN 

c= 6 
20 +n) (9) 


where k is the bulk modulus, 8 the cubical compressibility, u the rigidity, E 
Young’s modulus, ¢ Poisson’s ratio, and p the density. 


* This point is being investigated further. 
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The values obtained for these constants are: 


w=16.3+0.4 X10" dynes/cm? 
A=33+1X10" dynes/cm? 
k=44+1X10" dynes/cm? 

8 =2.28+0.05 X10-” cm?/dynes 
¢ =0.333 +0.005 
E=43+1X10" dynes/cm’. 


VI. Rockport REsutts 


Specimen records from Rockport are reproduced in Fig. 3 and the time- 
distance graph is shown in Fig. 5. 
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Fig. 4. Specimen Rockport records, vertical component. 


The velocity for longitudinal waves obtained from these data was 16,670 
+40 ft./sec. or 5.08 + 0.01 km/sec. 


"er 
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VII. Westerty REsvuLts 


The time-distance graph for Westerly is shown in Fig. 5. The velocity 
for longitudinal waves obtained from these data was 16,400 + 120 ft./sec. or 
5.00 + 0.04 km/sec. 
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Fig. 5. 


VIIL. AVERAGE LONGITUDINAL VELOCITY 


The weighted average of these three longitudinal velocities is 16,530 
+90 ft.//sec. or 5.04 +0.03 km/sec. 

The close agreement in velocity found here for three well-separated and 
different types of granite indicates that minor fluctuations in composition 
do not seriously affect the velocity. 


IX. ANALYSES 


Chemical 


The following analyses of specimens from the regions investigated were 
made by Mr. Forest A. Gonyer, chemist of the Harvard Department of 
Mineralogy. 


Mineralogical 


The following mineralogical analyses were made by Mr. Chalmer J.Roy, 
Assistant in the Department of Geology, who was a member of the party 
throughout the investigation. 

The Westerly analyses were made from thin sections. The Rockport and 
Quincy determinations were made by gravity separation methods from initial 
samples of 1000 grams. 
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TABLE I. 
Westerly Quincy Rockport 
Fine grain Coarse grain 
SiO. 74.00 70.30 72.86 74.37 
TiO, 0.10 0.31 0.18 0.18 
Al,O; 13.86 15.10 11.74 12.88 
FeO; 1.18 1.58 2.99 0.94 
FeO 0.71 0.99 1.56 1.56 
MnO none 0.04 0.03 0.03 
MgO 0.11 0.50 0.19 0.23 
CaO 0.75 L Ze 0.14 0.52 
NavO 2.88 3.69 4.36 4.11 
KO 6.20 5.58 5.28 4.92 
H.O — 0.09 none none none 
H.O4 0.31 0.40 0.41 0.31 
CO, none none none none 
P.O; none none none none 
SO none none none none 
100.19 99.74 99.74 100.05 
Taste IT. 
Westerly Quincy Rockport 
Fine grain Coarse grain 
Feldspar 45% 50% 60.10% 59.48% 
(orthoclase) (orthoclase) (mostly (mostly 
Perthite) Perthite) 
Quartz 35% 35% 28.88% 35.14% 
Others 20% 15% 11.02% 5.58% 
(muscovite 15 (muscovite 3 (riebeckite (hornblende) 
biotite 2 biotite 7 aegerite 
magnetite 3) magnetite 3 fluorite 
zircon and titanite) 
apatite 2) 


quartz and feld- quartz and feld- 


spar crystals spar crystals 
range from range from 
1/60 mm to 1g mm to 4 

2 mm in di- mm in di- 
ameter, aver- ameter, aver- 
age } mm ca. age 15 mm ca. 


X. Discussion 
La Courtine 


At La Courtine, Maurain, Eblé, and Labrouste® registered seismic waves 
generated by explosions and transmitted through a granitic terrane. 

They recorded longitudinal waves at distances of 5.6 km, 7.7 km, and 
13.9 km. When distances of this magnitude are used, a good velocity de- 
termination can be made with less accurate timing than is necessary with 
shorter distances such as were used in the present work. This advantage, 
however, is offset by the introduction of uncertainty as to the uniformity of 
terrane throughout the range. At La Courtine the presence of both gneiss and 
granite in the path of the waves lessened the certainty with which the veloci- 
ties obtained may be assigned to granite. 
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The method of timing the blasts introduced further uncertainty. The 
instant of the blast, because of the failure of an electric circuit designed to 
register it on the records, was determined visually by an observer equipped 
with a chronometer. An allowance of 0.2 sec. was made for his time of reaction. 

The mean longitudinal velocity obtained from five records in the La 
Courtine experiments was 5.524 km/sec. The five values, however, varied 
so much among themselves that they were segregated into two groups, giving 
averages of 5.905 km/sec. and 5.270 km/sec. respectively. The former was 
assigned tentatively to granite, and the latter to gneiss. An examination of 
their geologic map shows that this assignment is not valid because both 
granite and gneiss were traversed in every case. 

In spite of the use of “un amortissement trés faible”, the La Courtine 
investigators report, in addition to longitudinal waves, “autres, environ deux 
fois moins rapide, d’amplitude plus grande, se manifestent sur les trois com- 
posantes; nous les désignerons par le symbole ‘L’ appliqué de maniére 
générale en sismologie aux ondes de plus grande amplitude (longues ondes). 
De plus, a La Courtine, station la plus rapprochée des explosions, ont été 
enregistrés deux fois trés faiblement sur la composante horizontale trans- 
versale, et une fois en méme temps sur la composante verticale, des mouve- 
ments trés faibles oui paraissent ainsi correspondre a des ondes a vibrations 
transversale de vitesse intermediaire entre celles dex deux categories pré- 
cedentes; nous les désignerons par ‘S’ (secondes ondes).” The use of very 
feeble damping must, of course, have made the identification of these later 
phases somewhat uncertain. In the light of experience with transverse waves 
from explosions, which has been accumulated since the completion of this 
pioneering experiment, it seems clear that the waves designated as “L” 
were actually transverse. (See, for example, refs. 3, 7). The velocity of these 
waves, as determined from the mean of eight observations, was 2800 m/sec., 
making the ratio of the velocities of longitudinal and transverse waves 1.97. 
This is practically the same as the value of 2.00 obtained in the present in- 
vestigation. 

Thus, the velocities obtained by Maurain, Eblé, and Labrouste for both 
transverse and longitudinal waves, 2.800 km/sec. and 5.524 km/sec., re- 
spectively, are about 12 percent higher than the values for Quincy granite. 
Since both the method of timing the blasts and the presence of large areas 
of gneiss in the La Courtine region could introduce systematic errors, this 
divergence is not surprising. 


Comparison of laboratory and field results 


There is a fundamental problem involved in a consideration of the utility 
of comparing small-scale laboratory determinations of rock constants with 
field results. For many years it has been customary to use constants de- 
termined in the laboratory for computing elastic-wave velocities characteris- 
tic of certain rocks. On the basis of such comparisons deductions have been 
made concerning rock types represented by velocities observed in studies of 
near earthquakes. «t *” When the magnitude of the laboratory specimens in 
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past experiments is compared with the volume and field relationships of 
materials of which they are called representatives, however, this appears to 
be a very curious procedure indeed. 

A case in point is the investigation of Adams and Coker,! with particular 
reference to the granites which they measured, since those results bear di- 
rectly on the present work. Their measurements were made on granites from 
six widely separated localities. With the exception of Quincy, from which 
there were two, a single rock sample was used from each place, and from one 
to four test specimens were cut from each sample. The test pieces were about 
one inch in diameter and 3 inches long. The report states that “the rocks in 
all cases were air dry, having been allowed to remain in the laboratory for 
several weeks after they had been cut, before measurements were made”. 
Examination of individual measurements reported shows that constants 
vary from specimen to specimen, and in different directions in given speci- 
mens, up to as much as 50 percent. The lateral extension curve for one Quincy 
specimen has a definite kink, and Poisson’s Ratio at 2000 pounds of pressure 
is nearly 40 percent less than the average given for the specimen under pres- 
sures up to 9000 pounds. All of this raises a serious question as to the sig- 
nificance of the average results as representative of the Quincy granite in 
place. 

Laboratory determinations of the constants of a sufficiently large number 
of samples from a given region would, of course, give a much more repre- 
sentative figure. It is difficult to see, however, how past determinations re- 
ported in the literature can be said to represent general averages in any sense, 
entirely aside from any question as to methods and accuracy of measure- 
ments. Their acceptance as a basis for deductions from studies of near earth- 
quakes seems to imply an illusion that granite bodies, whether shallow or 
deep, are homogeneous masses of uniformly solid, flawless rock. Since they 
are anything but that, it should not seem strange if actual velocities of waves 
propagated through several thousand feet of granite in the field differed from 
those computed from the constants of isolated 1 in. X3 in. specimens of the 
same granite. 

Realizing this, we have nevertheless included the following comparisons 
on the possibility that they may record relationships between laboratory and 
field results which are interesting even if not particularly significant. 


Relation between compressibilty of granite and pressure 


Adams and Williamson? have determined the compressibility of several 
rocks and minerals by a static method at pressures of 2,000 and 10,000 mega- 
bars. In Fig. 6a graph showing the variation of compressibility with pres- 
sure is reproduced from their report. The dotted line and marked point 
have been added. It can be seen that according to this graph the statically 
determined compressibility of granite changes very rapidly at pressures under 
2000 megabars. Independent measurements of the compressibility of quartz 
and the feldspars which are the chief constituents of granite fail to show any 
such anomalous increase at pressures less than 2000 megabars, so the re- 












170 L. D. LEET AND W. M. EWING 


ported behavior of granite stands out as an exception to the law demon- 
strated by Adams and Williamson “that at moderately high pressures the 
compressibility of a fresh holocrystalline rock is an additive function of the 
compressibility of its minerals”, and resulted in their modifying it to apply 
only “provided the pressure is not too low”. The following quotations sum- 
marize the evidence which Adams and Williamson had for this change and 
show the difficulty which they encountered in explaining it. 

In the discussion on page 523 of their report, accompanying the graph 
shown in Fig. 6 here, is the following statement: “The results of Adams and 
Coker for six different granites provide the basis for the estimate of the range 
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Fig. 6. Compressibility as a function of pressure, according to Adams and 
Williamson. The circle and dotted line have been added to their graph. 


of variation for granites, at low pressures. The corresponding variation for 
the gabbros at low pressures is admittedly a very rough estimate representing 
mainly the opinion of the authors.” 

Thus, the only evidence for the anomalous behavior of the compressi- 
bility of rocks at low pressures seems to be the results of Adams and Coker 
on six granites. Even a futile attempt to find a consistent explanation for the 
anomaly apparently did not lead Adams and Williamson to doubt the validity 
of the earlier results as average figures for low pressures. Porosity was looked 
upon as a possible cause, but after considering it at some length, they con- 
cluded: “None of the rocks has enough porosity to explain, on the basis of 
detached voids or channels, more than a few percent of increase of com- 
pressibility. ... The granites, moreover, although containing less than one 
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percent of pore space, show an abnormally high compressibility at very low 
pressures; the value at a pressure of only a few hundred megabars being 
nearly double that at 1000 megabars.” 

After suggesting two more possible explanations and citing examples both 
contradictory and confirmatory, they proceed to a fourth: “Very often the 
abnormally high compressibility at low pressures is associated with coarse- 
ness of grain—for example, the granites and marbles, which are compara- 
tively coarse-grained, show a high initial compressibility. The Sudbury 
diabase, however, is moderately coarse-grained and does not show much of 
this effect, while the serpentine, which is very fine-grained, shows a con- 
siderable decrease in compressibility and this tendency to decrease persists 
even at high pressures”. 

In this same connection, it is interesting to note that Adams and Coker! 
had investigated the effects of grain size, including a granite in the trial runs, 
and concluded: “It will thus be seen that there is no correspondence between 
the coarseness of grain and the magnitude of the variations in the readings 
obtained.” 

This situation, which clearly proved rather difficult, and hinges essentially 
on the high values for the compressibility of granite obtained by Adams and 
Coker, is left by Adams and Williamson as follows: “About all that can be 
said concerning the way in which the compressibility of rocks changes with 
pressure is: For pressures above 2000 megabars the compressibility does not 
change very much; at pressures of only a few hundred megabars the com- 
pressibility is likely to be notably higher; and that the change of compressi- 
bility at low pressures, while connected with the admixture of minerals of 
different compressibility and with a looseness of structure existing in some 
coarse-grained rocks, cannot with certainty be predicted in advance.” 

The compressibility of Quincy granite, 2.2910~-" cm?/dynes, reported 
earlier in this paper, has an important bearing on the question under dis- 
cussion. The issue is somewhat confused by the fact that some writers con- 
tend that high-stress, statically determined compressibilities may be con- 
siderably different from those determined dynamically, which involve low 
stresses. If, for the moment, we assume the equivalence of statically and 
dynamically determined values, it is at once apparent that the compressi- 
bility of Quincy granite as plotted in Fig. 6 agrees with the high-pressure 
values of Adams and Williamson and eliminates the anomalous change at 
low pressures advocated by them. If, on the other hand, we deny this equiva- 
lance, it seems that it is almost inescapable to assume a variation in dynami- 
cally determined compressibility similar to that shown by Adams and Wil- 
liamson. The dotted curve in Fig. 6 represents this assumed variation. Ac- 
cording to this curve, the seismically effective compressibility at 2000 mega- 
bars would be roughly one half that observed at Quincy for pressures of a 
few megabars. Barring compensating changes in density and rigidity, such 
a decrease in compressibility would give a velocity of approximately 7.5 
km/sec. for longitudinal waves in granite at 2000 megabars pressure. This is 
far in excess of the figure 5.6 km/sec. observed by Jeffreys® and others as a 
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maximum velocity in the “granitic layer” to depths where the pressure is 
greater than 2000 megabars. We are accordingly forced to conclude that 
either the “granitic layer” is not granitic, or granite does not exhibit such 
rapid changes in compressibility at pressures below 2000 megabars. It ap- 
pears that, whether or not we consider statically and dynamically determined 
compressibilities equivalent, the observed compressibility of Quincy granite 
is not comparable with the low-pressure part of the Adams and Williamson 
curve. 

The good agreement between the compressibility measured at Quincy by 
a dynamic method and that measured statically by Adams and Williamson 
suggests that the marked discrepancy between static and dynamic values 
which has been reported*» does not exist for granite. It should be noted that 
though the only quantities which they measured were the compressibility 
and density, on page 520 of their paper they have computed the velocities 
of seismic waves in the rocks they studied. The values which they show for 
the velocities depend upon the value 0.27 which was assumed for Poisson’s 
ratio. If the value 0.333, which was determined at Quincy, were substituted 
in their computations, the resulting velocities would be in agreement with 
those determined in the present investigation. 
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ABSTRACT 

In the course of explorations of subsurface geology by the seismograph the authors 
have frequently observed the pronounced effect of stratification on the velocity of 
seismic waves in shales, and this effect has often been utilized in practical seismogra- 
phy. Recently an opportunity was afforded for securing additional quantitative data 
on the velocity normal to and parallel to the bedding planes. The paper points out 
that the velocity parallel to the planes of stratification is, in some instances, as much 
as fifty percent higher than the velocity in a direction normal to the bedding planes. 
It is shown also that inclined stratified beds exhibit a higher apparent point-to-point 
velocity when sound travels in an up-dip direction than when traveling down-dip. The 
paper describes a procedure whereby this effect may be utilized for determining the 
direction and approximate magnitude of the dip in such stratified deposits. The 
method has proved to be of considerable practical importance where the stratified 
formations are obscured by overlying deposits. 


GENERAL OUTLINE 


NE of the earliest applications of the seismograph was its use in profiling 
salt domes, limestone horizons, or other high velocity rock formations 
buried beneath sands, clays, or shales of lower velocity. In this method of 
profiling, the portion of the sound wave utilized travels substantially along 
the path shown diagrammatically in Fig. 1. As is well understood, the time 
interval between the arrival of the wave at two stations some distance apart 
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Fig. 1. Characteristic geology for diffraction profiling. 


along the line of propagation of the wave will be dependent on whether the 
wave is being propagated upslope or downslope. By measuring this velocity 
asymmetry the direction and magnitude of the slope of the high velocity bed 
can be determined. 

Commonly spoken of as diffraction or refraction shooting, this method 
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of profiling found wide application in many parts of the world, and still 
possesses much merit as a means of subsurface investigation in many local- 
ities. Its technique is now the most elementary part of the stock-in-trade of 
the oil seismologist. The procedure for obtaining a complete cross-sectional 
profile usually consists in shooting a line of detector stations from both direc- 
tions. Generally the overall time of travel of the sound wave is recorded from 
each shot to its appropriate stations, and complementary time-distance 
curves are plotted for the series of shots from each direction. 

Whatever may be the subsequent steps in making use of the primary time- 
distance data, all the calculations must necessarily be based on the deviation, 
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Fig. 2. Chart for showing relationships between apparent velocities when shooting upslope 
and downslope. Y =angle of slope; V; taken as 2500 m/sec.; V2 taken as 4500 m/sec. 


or asymmetry, of the apparent velocity relative to the true velocity of the 
rock which is being profiled. This asymmetry of the upslope and downslope 
values is best shown by comparing the station-to-station velocities derived 
by dividing the horizontal distance from one station to the next by the 
interval time, At. Fig. 2 shows the striking asymmetry which results from the 
slope effect in an ordinary case, while in salt dome work it is well-known that 
infinite apparent velocities are not uncommon. 

The above described method, as is well understood, gives the slope of the 
contact surface between the high velocity bed and the overlying low velocity 
deposits. In many cases this contact surface, where it has not been subject 
to irregular erosion, will conform to the structural outlines of the subsurface 
geology. In such cases the simple method above outlined is adequate to give 
the desired information on the subsurface geology. There are many cases, 
however, where this situation does not exist, and in such cases the method 
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fails to reveal the subsurface picture. Fig. 3 shows a typical case of this 
geology. Here the shale, 1, in which there exist more or less definite bedding 
planes, has been folded into a definite structural relief. Subsequent to this 
folding, extensive erosion took place leaving a surface shown by the line 2-3, 
and on this surface later deposits, 4, were laid down. In general, under these 
conditions, the shale, 1, practically always has a higher velocity than the 
recent deposits, 4. However, if we attempt to apply the simple method above 
described, we shall derive not the structural picture but merely a poor ap- 
proximation to the erosional surface between the shale and the recent de- 
posits. To further increase the difficulty there is often no definite reflecting 
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Fig. 3. Anticline hidden under unconformable beds. 


surface at greater depth which can be used to define the structure by reflec- 
tion shooting. Under such conditions the true structural picture cannot be 
revealed by any method of subsurface exploration heretofore published. These 
conditions are very frequently encountered and the purpose of this paper is 
to describe a method of exploration that is often capable of giving a correct 
structural picture where this situation exists. This method has been called 
“dip shooting” to distinguish it from the slope shooting heretofore used by 
geophysicists. 


ORIGIN OF Dip SHOOTING 


The method of profiling stratified formations which this paper outlines 
does not supercede any other method, but is simply a special development 
which is useful and applicable in certain cases which cannot be handled in 
other ways. The method is based upon an asymmetrical velocity effect which 
is analogous to the slope asymmetry which forms the foundation of refraction 
profiling. This “dip effect”, as it may be called, has been observed and util- 
ized in places as widely separated, geographically and geologically, as Canada 
and Venezuela. So far as our observations show, the phenomenon is common 
to most stratified rocks and even to some which are not ordinarily classed as 
stratified. 

Although the dip effect is analogous to the asymmetry produced by slope, 
it may arise from a different cause. The velocity asymmetry associated with 
stratification is often due to the effect of the bedding planes, dipping with 
respect to the surface, and not to a sloping subsurface. 

On this fact depend two important corollaries: First, the method neces- 
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sarily gives a structural rather than a topographic profile: for whereas the 
refraction method vields a structural profile under certain conditions, it is 
obvious that in the case of an eroded and buried land surface, for example, 
it merely records the old topography. Second, the method does not presup- 
pose the existence of any cover over the formation being investigated. In 
this case investigation with the seismograph sounds rather paradoxical, but 
the condition, nevertheless, is a real one and is frequently met with in 
tropical work. For example, shale areas which have been folded and sub- 
sequently leveled may be covered with just sufficient surface soil to make pit- 
ting difficult and expensive while heavy vegetation obscures all the outcrops. 
Here the seismograph, using the dip method, can profile the structure or 
work out the contours with almost as much facility as surface geology dis- 
plays in exposed regions. On the other hand, the presence of considerable 
recent overburden does not interfere with the use of the method. 

A third point which may be noted here as bearing on the theory of inter- 
pretation of the records, which we shall refer to later, is that the refraction 
system of profiling does not concern itself with the sound wave as it travels 
through the rock being profiled, but only with its emergence into the overly- 
ing formation. On the other hand, the dip method must take into account the 
phenomena attending the travel of the wave at a considerable depth within 
the body of the stratified rock itself. 


EVIDENCE OF VELOCITY ASYMMETRY 

The idea that stratification might influence velocity in a way to be of 
practical value was first suggested to the authors by J. E. Brantly. It re- 
quired considerable field worjgand the accumulation of extensive data, how- 
ever, before it was possible to sift the velocity variations due to stratification 
from the variations due to several other causes. Some of these other causes 
will at once occur to anyone who has used the seismograph. Among them we 
might list such factors as varying velocity in the overburden, dip in the over- 
burden compounding with the dip effect from the subsurface, differences of 
surface elevation, slope of the subsurface as distinguished from its dip in 
cases where the subsurface was eroded before being buried, ete. 

After making due allowance for such contributory causes as the foregoing, 
however, it was found that the sound wave travels faster parallel to the bed- 
ding planes than it does perpendicular to them. 

There are two lines of evidence to show that this is the case. The first is 
direct velocity measurement both parallel and normal to the bedding. 
Measurements of this kind were carried out recently on some exposures of 
the Lorraine Shales in the Province of Quebec, Canada. These shales have 
been highly folded along the Lake Champlain-Quebec Fault, and in a number 
of places are to be found standing on edge. This afforded an opportunity to 
measure the velocity perpendicular to the laminations of the shale and com- 
pare it with the velocity of the same formation in a horizontal position 
nearby. The following are some of the velocity readings. 
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TABLE I. Showing velocity of sound parallel and normal to bedding. 

















Velocity parallel to bedding Velocity normal to bedding 
planes. planes. Ratio 
(Meters per sec.) (Meters per sec.) 
4450 3205 1.39 
4655 3310 1.41 








From the above figures it will be seen that the velocity parallel to the bedding 
is in the neighborhood of 40 percent higher than it is in a direction at right 
angles. 

The second type of evidence which shows the existence of this asymmetri- 
cal effect is the difference in apparent interval velocity found in shooting 
up and down the dip of stratified rocks which have been tilted out of the 
horizontal position. Before presenting the field observations which illustrate 
this aspect of the problem, let us consider what theoretical explanation can be 
advanced to account for velocity variations as produced by dip. 

It is admittedly very difficult to make an exact analysis of wave paths in 
a non-homogeneous medium like shale, and we do not attempt to offer any 
comprehensive solution here. Even a preliminary study of the problem, how- 
ever, suggests that in addition to the refraction produced by more or less 
alternating layers of different velocity, there is also refraction due to increase 
of velocity with depth in most cases. Furthermore, it seems logical to suppose 
that considerable absorption of the wave energy would occur when the wave 
path cuts across the bedding, and that this absorption would tend to dis- 
appear when the path parallels the stratification. 

Consideration of these and similar factors affecting the travel of the sound 
wave leads to the conclusion that we can makéa useful approach to the prob- 
lem by studying some of the simpler typical subsurface conditions such as 
are often met with in the field. From these simple cases it should be possible 
to extend the analysis to more complicated conditions with the aid of further 
experimental work. It is to be hoped, therefore, that additional data bearing 
on the subject will be forthcoming from localities where the dips are more 
accurately known than in many areas which the authors have investigated. 

To return, however, to those instances where it can be readily seen that 
dipping strata would give rise to asymmetrical velocities, the three clearest 
cases, perhaps, are the following: 


CASE 1 


The first example in which dip will obviously give rise to velocity asym- 
metry is illustrated diagrammatically in Fig. 4. Almost any stratified forma- 
tion has occasional hard layers which give a somewhat higher velocity than 
the average. These layers may be of the nature of “markers” which are readily 
distinguished by some peculiarity, or they may be almost unrecognizable to 
the eye while still possessing abnormal velocity characteristics. Such layers 
act as conductors for the sound wave because they offer the path of least 
time. 
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In the diagram, Fig. 4, layers of the above-the-average velocity are repre- 
sented at a-a. In shooting a line of stations in an up-dip direction, the time 
intervals between stations 1-2 and 2-3 will be short, while between stations 
3-4 and 4-5 the intervals will be longer, because these stations are beyond the 
zone of influence of the fast layer “a”. However, on advancing farther down 
the line, other laminations of higher-than-average velocity will be en- 
countered and the intervals will again be shortened. Conversely, when shoot- 
ing down-dip, many of the intervals will be lengthened by these same layers. 

Now it is expectable that the flatter the dip (until it becomes nearly 
horizontal), the more continuous will be the effect on the time intervals. 
Therefore, moderate dips, such as are usually found in practice, will show this 
type of velocity asymmetry to the best advantage, and those formations 
which exhibit the least uniform stratification will generally be the ones whose 
dip it will be easiest to determine. 
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Fig. 4. 


It may be objected that the foregoing analysis is based on the slope- 
effect and not on dip. To a certain extent this is true. In extreme cases, like 
that shown in Fig. 4 where the high velocity layers are unusually thick and 
clearly defined, the interpretation can be based purely on the theory of 
refraction shooting. But in a majority of cases the layers of higher velocity 
are relatively thin, rather closely spaced, and seldom recognizable in the 
outcrops. Thus it seems to be the aggregate effect of many sloping layers 
which generally produces asymmetrical velocities such as those tabulated in 
the latter part of this paper. Consequently, there is no choice but to make the 
interpretation in terms of dip rather than in terms of slope. 


CASE 2 


In Fig. 5a the wave is shown traveling in an up-dip direction. Exaggerat- 
ing conditions somewhat for the sake of clearness, we may assume that the 
wave, taking the path of least time, descends fairly deeply into the shale and 
is steadily refracted on a curved path until the critical angle is reached, when 
it follows the bedding and finally emerges at the surface giving the time 
interval Af,. In shooting down dip, the principle of reversibility of wave paths 
tells us that, between identical surface points, the path will be the same from 
either direction and there can be no difference in the recorded times. Simi- 
larly, if the subsurface dip is everywhere uniform, equal shooting distances 
will give identical paths even though the surface points are not the same. 
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Krom this it follows that a perfectly uniform dip, in this sense, will not result 
in any velocity asymmetry. 

On the other hand, if the dip is non-uniform as illustrated in Fig. 5 and as 
is most commonly the case under actual conditions, especially close to an 
anticline, the wave paths will not be the same provided the surface points 
are shifted. Thus, in Fig. 5b, the shooting distance, L, is the same as in 5a 
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but the setup is shifted slightly so as to intercept different conditions of dip. 
(“onsequently, the wave paths are not the same for the two cases, and both the 
overall times and the station-to-station time intervals will be different. In 
other words, Af, will not, in general, have the same value as at At, and the 
apparent velocities will show an asymmetrical effect. 
CASE 3 

A third type of velocity asymmetry is often due to a relationship which 
frequently exists between the hardness and the age of rocks and shales. In 
shooting up-dip (except in the case of overturned folds), the wave is traveling 
from younger to older beds. Now in many cases it has been found that the 
hardness of rocks increases with their age and the velocity closely corresponds 
to the hardness. Consequently, in shooting from a given shot point to a de- 
tector setup located in an up-dip direction, the tendency is for the velocity 
to increase, and the converse is true in shooting down-dip. This tendency, 
while by no means invariable, must nevertheless be taken into consideration 
as a causative factor. Whenever this effect exists it gives a true indication of 
the direction of the dip of the bedding planes. 
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The reader’s attention has already been called to the complex nature of 
the problem under discussion. It is well to emphasize again, however, that 
many factors affecting the velocity of sound in stratified rocks have not yet 
been mentioned. Indeed, many factors doubtless exist which have not yet 
been recognized. Therefore, there is still room for much experimental in- 
vestigation along these lines. At the present time, lacking a theoretical ex- 
planation which covers all phases of the phenomenon, the existence of the dip- 
effect and the validity of its use as a working method are best evidenced by 
numerous field observations; a few of which follow. 


EXAMPLES OF THE Dip Errect 


TABLE I]. Material: Tertiary shale in Venezuela. Shooting distance: 1500 meters. Station spacing: 
250 meters. 





NW Shots SE Shots 











Stations Velocities Stations Velocities 
442-411 2450 106-107 2660 
111-110 2360 107-108 2720 
110-109 1980 108-109 3120 
109-108 2080 109-110 2600 
108-107 2320 110-111 2500 
107-106 2400 111-112 2720 
106-105 2400 112-113 2200 
105-104 2600 113-114 2180 
Average 2324 2588 








The velocity values shown in Table II were recorded in the course of 
surveying a line which was laid out approximately at right angles to the 
general strike of the area. The total distance covered by the stations is two 
kilometers, and the stations have been so chosen as to represent the same 
portion of the subsurface in the northwest and southeast shot groups, after 
making allowance for the probable amount of throwback of the reference 
point. The shale, in this case, was covered by clay, but the evidence of various 
nearby test wells pointed to the cover being fairly uniform in thickness. Con- 
sequently, the difference in the velocity averages indicates a regional dip to the 
southeast. The variations of the individual velocity values show the effect 
of local folding. 

We have already mentioned the fact that some relatively unstratified 
formations show asymmetry due to dip. The following data were gathered in 
areas of indistinctly bedded sands and clays in Western Venezuela. 

The values in Table III cover a distance of ten kilometers with an allow- 
ance for throwback as in the case of Table II. The difference between the 
average velocities is small, but in view of the circumstances a larger differ- 
ence would hardly be expected. These shots were taken at a distance of about 
twenty kilometers out from the foothills of the Perija mountains. The dip 
represented by the averages is simply the normal depositional dip of water 
borne sediments laid down along the mountain front. The dip is toward the 





182 B. McCOLLUM AND F. A. SNELL 


Tase III. Material: Tertiary sandstones, sands and clays. Shooting distance: 6000 meters. 
Station spacing: 500 meters. 








E Shots W Shots 
Stations Velocities Stations Velocities 
165-163 3290 127-129 3420 
163-161 2980 129-131 2980 
161-159 3120 131-133 3090 
159-157 3120 133-135 3050 
157-155 3050 135-137 2970 
155-153 3420 137-139 2990 
153-151 3380 139-141 2940 
151-149 3160 141-143 2910 
149-147 3050 143-145 3010 
147-145 3900 145-147 2600 
145-143 2940 147-149 3730 
143-141 3210 149-151 2870 
141-139 3290 151-153 3470 
139-137 3180 153-155 2910 
137-135 3080 155-157 3120 
135-133 3090 157-159 2980 
133-131 2930 159-161 3290 
131-129 3080 161-163 3160 
129-127 3080 163-165 3200 
127-125 3290 165-167 3330 
Average 3182 3101 








east, that is, away from the mountains, and would doubtless be steeper if 
investigated nearer the foothills. To illustrate this latter point we may 
average the western half of the group of velocities, Stations 125 to 147, with 
these results: East shot average, 3117 meters per second; West shot average, 
2996 meters per second; showing a greater difference than for the whole 
group. 


TABLE IV. Material: “Mottled clay”, Venezuela. Shooting distance: 500 meters. Station spacing: 
500 meters. 


————————————— ———— 











E Shots W Shots 
Apparent velocities Apparent velocities 
2840 2480 
2870 2480 
2340 2500 
2530 2250 
2810 2290 
2400 2420 
2780 2170 
2600 2140 
Average 2646 2341 








Table IV shows the results of some scattered shots on two lines at right 
angles to the Perija Mountain front. The material is what is known as the 
“Mottled clay” of the Lake Maracaibo region. It is a compact clay with very 
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little admixture of sand and has no recognizable bedding. In the locality 
where the above data originated, however, there were good reasons to suppose 
that the regional dip was eastward, in a direction away from the hills. Drilling 
subsequently indicated that such was the case for the underlying formations, 
but the seismograph results appeared to be the only definite evidence that the 
same dip existed in the mottled clay cover. Shots at 1000 meters distance 
confirmed the figures obtained from the 500 meter shots, as the following 
tabulation will show. 


TaBLE V. Material: “ Mottled clay”, Venezuela. Shooting distance: 1000 meters. Station spacing: 
500 meters. 











E Shots W Shots 
Apparent velocities Apparent velocities 

2320 2360 

2660 2550 

2780 2550 

2690 2580 

2940 2580 

2800 2580 

2230 

2170 

2320 

Average 2698 2436 














Tests ror KNown SrRUCTURES 
Perhaps the best confirmation of any geophysical theory is its test on 
known structures. We have carried out several tests of this sort with very 


gratifying results. One of the earliest, and at the same time, one of the most 
comprehensive trials of the dip method of shooting was made on a group of 
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Fig. 6. Dip-velocity curves across structure miocene shale—Venezuela. 


shale structures, arranged in echelon, near Maracaibo, Venezuela. These 
folds, outcropping through the mottled clay cover, had been worked out by 
surface geology, and their axes fairly well located. The seismograph was 
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employed to check the location of the recognized folds and to search for others 
which were not exposed on the surface. 

The results tallied with the geological picture in a manner which could 
hardly be questioned. It is not possible, here, to present all the data which 
were gathered, but one example may not be out of place to illustrate this 
special application of the method to exposed or very thinly covered folds 
which have been planed by erosion before being buried. Fig. 6 shows the 
apparent velocities given by shots from opposite directions plotted in the form 
of velocity curves. The line of survey crossed one of the known shale struc- 
tures already mentioned. The cross-section of the fold is shown diagram- 
matically at its proper location near Station 106. The northwest flank of the 
anticline is considerably steeper than the southeast flank. In shooting from 
the NW, note how the velocity rises to an abnormally high value due to the 
steep dip of that flank. A corresponding minimum value is given by the south- 
east shots, and the reverse occurs in both cases on the other side of the fold 
after crossing the crest. Doubtless a part of the velocity asymmetry is due to 
slope effect on the surface of the shale under the clay cover, but repeated ob- 
servations seem to demonstrate that this will not account for all of the 
velocity variation which is found in such cases. 


APPLICATION OF Dip SHOOTING 


The main essential in utilizing the dip effect is to shoot all stations from 
opposite directions. Generally the stations are set along a straight line which 
is laid out, as nearly as may be, at right angles to the strike of the area in 
order to get the maximum effect from the dip. Shooting distance, station 
spacing, etc., must be determined according to what is known of the surface 
and subsurface conditions, as with any other method. Often special prelim- 
inary shooting is necessary to reveal this information. 

It is possible that satisfactory results might be obtained from the use of 
overall velocities, derived from the total time of travel of the wave from shot 
to detector, but the authors have always preferred to employ station-to- 
station velocities because they not only localize the dips better, but also bring 
out the velocity asymmetry more sharply. 

In accordance with this practice, then, the interval velocities given by 
the shots from both directions may be tabulated and the complementary 
values compared. It is necessary to make allowance for the throw-back of the 
reference point as is customarily done in all profile work, although here the 
exact amount to allow for is indeterminate on account of the variable form 
of the wave path. The fact that the wave may penetrate to an appreciable 
depth in the shale, and that all parts of the path are generally curved, 
renders it impossible to do more than make an approximate calculation of the 
throw-back, but this will usually suffice for practical purposes. 

The velocity asymmetry observed on any setup may, of course, be due to 
two causes. First, the slope of the contact between the shales and the overly- 
ing formations, and second, the velocity asymmetry in the shales themselves. 
In order to determine and eliminate the former, it is usually best to shoot a 
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correction shot on each setup from a point symmetrically placed with respect 
to two detector stations, as, for example, at a point midway between them. 
If the spacing of the stations has been properly chosen with respect to the 
local geological conditions, which can always be accomplished, the required 
correction can easily be made. 

The most difficult phase of the work is to translate the velocity readings 
into terms of degrees of dip. There seems to be no theoretical basis for a 
rigorous mathematical treatment which would accomplish this. Consequently, 
the best solution is probably an empirical calibration of the velocity values 
on some structure where the dips are accurately known. In most cases, how- 
ever, it is not of great importance to know the absolute dip. For general 
profiling and reconnaissance the direction and relative steepness of the dips 
is all that is required. The use of relative values permits comparison of the 
size of different structures, shows which side of the fold is steeper, and, of 
course, locates the crest of the fold as accurately as if precise quantitative 
dips were obtained. 


SUMMARY 


A large amount of accumulated data similar to that herein presented 
demonstrates that the method described above provides a satisfactory means 
of outlining structure in stratified and other asymmetrical deposits where 
conditions are such that the usual methods of refraction and reflection cannot 
be applied. 
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ABSTRACT 

The equation of motion of a mechanical seismograph is —%=+2kj+p*y where 
x is the ground displacement and y the seismograph deflection. This equation may be 
solved for y when x is supposed known or for x when y has been observed as a function 
of the time. In this paper both of these ways of solving the equation are considered. 
The motion of the seismograph due to a train of waves starting at t=0 is considered 
and also the motion due to the arrival of a single wave. In each case seismographs with 
several periodic times and either undamped or critically damped are considered. 
Curves are given showing the motion of the ground and the calculated motion of the 
seismograph. The motion of the ground corresponding to several simple assumed seis- 
mograms is also worked out and shown by means of curves. The motion corresponding 
to a given seismogram depends greatly on the periodic time and damping of the seis- 
mograph., Finally the ground motion is deduced from two actual seismograms due to 
dynamite explosions. An integraph is described which enables the calculations to be 
done more quickly. 


N THE case of seismograms obtained with mechanical seismographs hav- 

ing a definite periodic time and damping coefficient the motion of the 
ground can be deduced from the seismograms without serious difficulty. This 
of course is well-known, but so far as the writer is aware few if any results of 
such calculations have been published, at any rate for seismograms obtained 
in geophysical prospecting.! 

The seismograph consists essentially of a heavy mass supported on 
springs, or in some other way, so that it is free to oscillate. The motion of the 
ground relative to this mass is magnified and recorded. 

Consider a mass m supported by a spiral spring from a stand resting on 
the ground. Let y denote the upward displacement of the mass m, relative 
to the stand, from its equilibrium position, so that y is equal to the decrease 
in the length of the spring due to the upward displacement of the mass relative 
to the stand. When y is zero the force exerted by the spring is equal and op- 
posite to the weight of the mass. 

Let x denote the upward displacement of the ground and stand due toa 
wave in the ground relative to axes supposed fixed. The equation of motion 
of the mass m is then 


mi + ¥) = — mp*y — 2mky. 


1H. Arnold, Zeits. f. physik. Erdkurde 10, 269-317, gives calculations of the ground mo- 
tion from several seismograms. 
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Here p” is the restoring force per unit mass exerted by the spring when its 
length is increased by unity and 2k is the damping force or viscous resistance, 
per unit mass, to the motion of the mass m when moving relative to the stand 
with unit velocity. Hence we have 


i+ 22kpt py = — 4. 


We see from this equation that if initially x=0, =0 and the mass is at 
rest in its equilibrium position, then when the ground begins to move we shall 
have y= — x because at the start y and j are zero so that j = —#. This merely 
means that at the start the mass m remains at rest relative to the fixed axes. 
But y will only be equal to —x during an interval small compared with the 
natural period of vibration of the mass mz. 

The above equation may be regarded as giving y when x is known or as 
giving x when y is known. In previous discussions of the theory of seismo- 
graphs, which the writer has seen, the motion of the ground or x has been as- 
sumed known and the resulting values of y have been calculated. 

The solution of the equation j+2ky+p*y = —¢= F(t) where F(t) denotes 
a function of the time ¢ and when y and ¥ are zero at t=0 is 


1 t 
_ —{f e~*(t—”) sin p(t — n)F(n)dn 
Bh Jo 


where 


m (p? _— Rk?) 1/2, 


By means of this equation it is easy to calculate the seismograph deflec- 
tions y as a function of ¢ for any simple values of F(¢). 

In homogeneous ground an explosion of dynamite would be expected to 
produce single waves, or short trains of waves, so that it is important to con- 
sider the motion of the seismograph due to single waves. To represent a single 
wave we may suppose that x =A(1—coswt) from t=0 to t=27/w and that 
x =0 when ¢<0 and when t>27/w. This gives the wave shown in Fig. 1. 





Fig. 1 Fig. 2. 


Before considering this case it will be convenient to consider the case when 
x =A(1—coswt) when t>0 and x=0 when ¢<0 which gives the series of 
waves shown in Fig. 2. 

The equation x = A(1—cosw#) gives ¢= Aw’ cos wt so that the expression 
for y becomes 


Aw* ¢! ; 
y= -— e~*() sin w(t — ) cos wydn. 
w Jo 
We shall consider two cases. One when the seismograph is undamped so that 
k =0 and one when the seismograph is critically damped so that »=0. 
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In the first case with k=0 we get 


Aw? 
y = —(cos wt — cos pt) 
w? — p? 





and in the second case with np =0 we get 








Aw? Aw*e—** 
y = — cos (wt + a) + - ———1}$k?(1 + kt) — w(1 — kt) } 
w? on 2 (w? oa k?)? 
where 
2wk 
tana = — — 
w? — k2 


We shall now consider three cases with the undamped seismograph (a) 
when the period of the seismograph is greater than that of the waves in the 
ground (b) when the period of the seismograph is equal to that of the waves 
and (c) when the period of the seismograph is less than that of the waves. 

In case (a) let w=27/12, p= 27/36 and A =1 so that 


y = — (8/9)(cos (304)° — cos (10#)°) 
x = 1 — cos (30/)°. 


Fig. 3 shows the curves given by these equations. We see that the movement 
of the ground sets up an oscillation of the seismograph so that the deflections 





Fig. 3. 


obtained are the resultants of the oscillations of the ground and of the seis- 
mograph. 





Fig. 4. 


Fig. 4 shows the curves obtained with w=27/12 and w=27/120 so that 
y = —cos(3t)°+cos(30t)° « = 1 —cos(302)°. When, as in this case, the period of 
the seismograph is ten times that of the waves, the undamped seismograph 
records the first wave fairly accurately. 

Now consider case (b) when the period of the seismograph is equal to that 
of the waves. If we put w=y-+e in the equation 
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Aw? 
y = ——(cos wt — cos pb) 
ee 
we get when € is very small y= —Auyt/2 sin wt. If A=1 andw=yu=27/12 this 
gives 
* t sin (301)° 
—y=—tsin (; 
al (304) 
and 
x = 1 — cos (30#)°. 


Fig. 5 shows the curves for y in this case. The amplitude of the first 
oscillation is only about one half that of the ground but the amplitude rapidly 
increases as the successive waves arrive. 








Fig. 5. 


We see from this that for the purpose of detecting the arrival of the begin- 
ning of a train of waves, it is not a good plan to use a seismograph having a 
free period equal to that of the waves. A seismograph with a longer period 
gives twice as great a deflection due to the first wave as one with a period 
equal to that of the waves. Of course at the start — y= x in all cases so that 
the initial sensibility is independent of the periods but it is advantageous to 
get a larger deflection due to the first wave since a large deflection is less 
likely to be overlooked or obliterated by the small oscillations due to local 
disturbances. 

Now consider the third case when the period of the seismograph is smaller 
than that of the waves. The equation 


Aw? , 
7° —— (cos wt — cos wt) 


pw? — w? 
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shows that making yu/w large will diminish the maximum values of y very 
greatly. For example if u/w=10 we get 


y = (A/99)(cos 10wt — cos wt) 


so that the maximum amplitude is only about one percent of that of the 
waves in the ground. In this case we get short oscillations due to the term 
cos 10wt superposed on the curve given by —(A/99) cos wt. Hence a smooth 
curve drawn across the small oscillations would give —(A/99) cos wt which 
resembles the curve given by x = A(1—coswt) although displaced downwards 
and with only one percent of the amplitude. This suggests that the results 
obtained with high frequency undamped seismographs can be improved by 
drawing a smooth curve through the small oscillations of frequency equal to 
that of the seismograph. 
Fig. 6 shows the curves given by 


x = 1 — cos (302)° 


and 


— y = (})(cos (90#)° — cos (308)°) 


which are those for the case A = 1, w=27/12 and w= 27/4 in which the period 
of the seismograph is one third that of the waves. 





Fig. 6. Fig. 7. 


Now consider the case of a critically damped seismograph. In this case if 
w= 27/12 so that x = 1—cos(30t)° and k= p= 27/36 so that the seismograph 
period is three times that of the waves then the equation for y gives 

9 Sat — 72 
y = — cos (30% + 36-87)° fe ne go POIS 
10 100 

Fig. 7 shows the curves given by these equations. Comparing Fig. 7 with 
Fig. 3 we see that the critical damping cuts out the oscillations of the seis- 
mograph, set up at the start, so that we do not get interference between two 





Fig. 8. 


vibrations of different frequencies. The damping however diminishes the 
interval at the start in which x and —y agree. Thus without damping in Fig. 
3 the two curves agree from ¢=0 to ¢=4 but with critical damping in Fig. 
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7 they only agree from ¢=0 to t=1. We see that to obtain the actual motion 
of the ground it is necessary to use a much longer period with a critically 
damped seismograph than with an undamped one. 

Fig. 8 shows the curves for x and —y when w= p= so that the periods of 
the waves and of the seismograph are equal. The equations for x and y in this 
case are 
x = 1 — cos (301)° 
y = — 3 sin (304)° + (at/12)e—**/6. 


The seismograph is less sensitive and x and —y do not agree as well as when 
the period of the seismograph is greater than that of the waves. 





loy 
Fig. 9. 


If the period of the seismograph is one third that of the waves and the 
damping critical then we may take 


x = 1 — cos (307)° 
so that w=7/6 and k=7/2. The equation for y is then 


1 Sart + 8 
y = — — cos (30 — 36-87)° + ———_ 2, 
10 100 
Fig. 9 shows the curves for x and —10y in this case. —10y is plotted instead 
of —y because the maximum value of y is only one twentieth that of the 
maximum value of x. 

Comparing Fig. 9 with Fig. 7 we see that the deflections of a critically 
damped seismograph with period three times that of the waves are somewhat 
similar in form but about ten times larger than those of a critically damped 
seismograph with period one third that of the waves. 

It appears that, while in all cases the initial deflections are equal and op- 
posite to the motion of the ground, the only case in which this remains true 
approximately for more than a fraction of a wave is that in which the period 
of the seismograph is large compared with that of the waves and the seis- 
mograph is undamped. 

The ratio of the maximum value of y to the maximum value of x is about 
w*/(w*—y*) or if T is the period of the waves and 7” that of the seismograph 
itis 7’*/(T’*—T*). This is equal very nearly to unity when 7’/T is large and 
is nearly equal to (7’/T)? when 7’/T is small. When T’ and T are nearly 
equal it becomes large but not for the first wave. 

Very short period seismographs are difficult to make sufficiently sensitive 
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because of the small value of the ratio of the maximum value of y to that of 
x. A critically damped seismograph with a period very small compared with 
that of the waves in the ground gives a deflection very nearly proportional 
to the acceleration of the ground. The equation 
1 t 
y=- i) e~*(9 sin w(t — n)F(n)dn 
ab Jo 


when » = 0 for critical damping becomes 


{ ek (t — 9) F(n)dn. 


If the period of the seismograph is very small so that & is large then the fac- 
tor e~*'-" is negligible except when 7 is nearly equal to ¢. If we suppose that 
F(n) varies slowly and so can be regarded as constant during the short inter- 
val when 7 is nearly equal to ¢ then 


t 
y = Fié) f e~*'—-(t — n)dn 
0 


or 
F(t) 
y = — (1 — o **(1 + Re)) 
kb? 
or when ¢ is not very small y= F(t)/k?. Thus in this case since F(t) = —# we 
have y= —i#/k’. 


If a very short period critically damped seismograph can be made suffi- 
ciently sensitive, it will indicate the acceleration of the ground with consi- 
derable accuracy. Whenever a sudden change in the acceleration of the ground 
occurs, the change in the seismograph deflection will be proportional to the 
change in the acceleration multiplied by the factor 1—e-*‘(1+¢) where t¢ is 
the time since the sudden change took place, which very soon becomes equal 
to unity when & is large. 

Now consider the deflection of the seismograph due to a single wave. Let 
x =A(1—cos wt) from t=0 to t=27/w and x =0 when ‘<0 and when ¢> 27/w. 
In this case y is given by 

Aw ' 
y= ——] e&*™ sin w(t — 9) cos wndy 
M 0 


from t=0 to {=27/w and by 


ef fu 
y= -— e~*'™ sin w(t — 9) cos wndn 
MB 0 


when / is greater than 27/w because the acceleration is zero when ¢>27/w. 
First suppose that the seismograph is undamped so that k=0 and then 
when />27/w we have . 
Aw?* 2r/w 


y= -— sin u(t — n) cos wndn 
MB 0 





alia 
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which gives 
2Aw* sour 
y = —— sin — sin p(t — w/w). 
w? — pb? w 
From t=0 to ¢=27/w the equation for y is 
Aw?* 
y = ——~(cos wt — cos ut) 
w? — p’ 
as for x = A(1—cos wt) when t>0. 
The curves for y and x in the three cases 
A=1, w=7/6, np=7/18 
A=1, w=7/6, p= 7/6 
A=1, w=7/6, p=2/2 


are shown in Fig. 10, Fig. 11, and Fig. 12. 





Fig. 10. 








Fig. 11. Fig. 12. 


In Fig. 10 the period of the seismograph is three times that of the wave 
and a large oscillation is started in the seismograph. In Fig. 11 the two periods 
are equal. In Fig. 12 the period of the seismograph is one third that of the ' 
wave and the wave ends with the seismograph in equilibrium so that y=0 
after the wave has passed. 
Now suppose that the seismograph is critically damped so that »=0 
and k=p. In this case the equation for y when ¢>27/w is of the form 
y=at+b(t—2r/w)e—*'**/) and the curves from t=0 to t=27/w, of course, 
are the same as with x = A(1—cos wt) when t>0. Figs. 13, 14 and 15 show the 
curves for —y and x in the three cases 
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A=1, w=7/6, p= 2/6 


A=1, w=7/6, p=7/2. 


We see that with critically damped seismographs the shape of the curves ob- 
tained depends less on the period of the seismograph than with undamped 
seismographs. When the period of the seismograph is much smaller than that 
of the wave the deflections are very small. 








lOy 


Fig. 14. Fig. 15. 


It appears that the only case in which the deflection of the seismograph is 
approximately equal and opposite to that of the ground during the passage 
of one or more waves is that of an undamped seismograph with period much 
larger than that of the waves. In this case, however, a single wave sets up an 
oscillation which continues indefinitely. 

The principal conclusion which can be drawn from the above results is 
that seismographs do not register the actual motion of the ground so that 
arguments based on wave forms or the apparent arrival of waves later than 
the first wave must be used with great caution. 

Critically damped seismographs are obviously better when it is desired to 
observe the arrival of several waves separated by appreciable intervals. 

The inverse problem of calculating the motion of the ground from the 
recorded deflections of the seismograph will now be considered. The equation 
—#=j+2ky+p’y on integrating gives —i=yt+2kyt+p?f ydt+C. We may 
suppose that y=0 at =0 so that —X%»=so+C. But at t=0 —xX9=%0 so that 
C=0. Also when tis very large we may suppose * = j = y =0 so that /o”ydi =0. 
Integrating again we get 


t t t \ 
—x=yt+ 2k f ydt + pf § f ydt bat + D. 
: ie F 


If x=y=0 at t=0 then D=0O. Let J’ydt =a, and Sh fiydt} dt =a2 so that 
—x=y+2ka,+ p*a,. When ¢ is large we may suppose that x=y=0 so that 
we must have 2ka,+ p*a.=0. But a, =0 when t= © so that a2=0 whent=o. 
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In the equation —x=y+2ka,+ p*az since p? =472°/ 7? where 7 is the natural 
period of the seismograph, when undamped, we see that the correction pa, 

which must be added to y to get the true motion of the ground is inversely as 

T? for all values of the time ¢. If 7 is very large this correction may be negligi- 

ble. In the same way the damping correction 2ka; is proportional to the damp- 

ing coefficient. Thus for an undamped seismograph of very long period —x=y 

so that the seismograph gives the ground motion exactly. Any increase of the 

damping and any reduction of the period 7 increases the corrections and so 

increases the difference between the seismograph deflections and the ground 

motion. If k= the seismograph is critically damped. 

To find a, it is convenient to make a tracing of the seismogram on graph 
paper and draw a curve for the area by counting the squares between the y 
curve and the zero line y =0. The zero line is usually not marked on the seis- 
mograms so that it is necessary to draw it. A small error in the position of the 
zero line makes a large error in a;. Thus if it is drawn at y = —} instead of at 
y=0 we get a, too large by 3/. Such errors usually occur but are easily allowed 
for when it is remembered that a,;=0 for large values of ¢. If a; does not fi- 
nally approach zero as ¢ increases then by slightly shifting the zero line it is 
made to do so. Very often it is found that the area curve fluctuates about a 
straight line through the origin finally approaching this line. In such cases it 
is clear that the zero line is too high or too low. The line through the origin 
may then be taken to be the zero line for the area. The integral a2 is then ob- 
tained from the curve giving the area in the same way. If it does not finally 
approach zero it must be made to do so by a suitable shift of its zero line. If 
an error Aa is made in the area a; at time ¢, then this produces an error 
Aa(t—t;) in the value of a2 at a time ¢ later than ¢;. Such errors are usually 
appreciable and cause the curve obtained for this integral to drift up and 
down. It is therefore usually necessary to draw a new zero line for this curve. 
This zero line must be a smooth curve such that the value of the integral 
oscillates more or less equally above and below it. 

Instead of endeavouring to correct such errors by shifting the zero lines, 
a process which depends a good deal on arbitrary judgment, one may carry 
out the calculation of x using the uncorrected curves for the area and the 
integral. When this is done the curve obtained for the motion of the ground 
usually appears to oscillate about a smooth curve which may drift very far 
from the zero line. This smooth curve may be drawn and taken to be the 
true zero line for the ground motion. These methods of correcting the calcu- 
lated ground motion are necessary unless the areas are very accurately de- 
termined. They obviously result in the removal from the ground motion of 
any long period oscillations having periods long compared with the period of 
the seismograph. The constant p is nearly equal to 27/T unless the damping 
is large, T being the periodic time of the seismograph. In getting the areas of 
the curves it is convenient to use a time scale such that p=1. 

Before considering results obtained with seismograms some simple the- 
oretical cases will be considered. Suppose y=0 when ¢<0 and y=a sin pt 
when ¢>0 and that k =0. In this case 
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a sin p+ op f | f sin pe ath a 
0 0 


= asin pt + ap { (1 — cos pt)dt 
0 


| 
I 


which gives 


z 
| 


so that 


— x = asin pt + pat — asin pt = pat. 


Thus in this case the ground moves with uniform velocity ap. The sudden 
starting of this motion at t=0 starts the seismograph oscillating freely with 
amplitude a but has no further effect on the motion. If a ground wave ar- 


y 





k=O, T=7/2 





Fig. 16. 


rives the period of which is long compared with that of the seismograph it 
will produce a nearly uniform motion of the ground for a short time which 
will start the seismograph oscillating with its natural period. 
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Suppose that y=0 when ¢<0O and y=a cos pt when ¢t>0 and k=0. In this 
case —x=a cos pt+ap?fo'{ fot cos ptdt}dt which gives —x=a. Thus in this 
case the ground suddenly moves a distance a at t=0 and then remains at rest. 











0 Te2n  k4/2 O K2  T=2n 
Fig. 17. Fig. 18. 


Figs. 16, 17, 18, show the curves for y, ai, a2, and —x in some simple the- 
oretical cases which need no explanation. The seismograph of course always 
records sudden movements correctly. Fig. 19 shows a seismogram of the 
waves due to an explosion of 1000 pounds of dynamite at a distance of 9870 
meters. The seismograph was of the vertical mechanical type. The calculated 
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Fig. 19. 


curve for the ground motion (x) is also shown. It appears that the explosion 
waves consist of a rapidly damped oscillation of period about 0.1 sec. Several 
such waves come in at different times and all have about the same shape and 
period. In this case there were no geological anomalies in the neighborhood, 
but the wave velocity increased more or less uniformily with the depth. The 
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waves therefore travelled along nearly circular paths and may have been re- 
flected from the ground surface as shown in Fig. 20. In this way the arrival of 
several waves at different times may be explained. Fig. 21 shows a seismo- 
gram of the waves from 1000 pounds of dynamite at 9590 meters. In this case 
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Fig. 21. 


there was a salt dome in the path of the waves. The curve for x the ground 
motion is also shown. The ground motion can be regarded as mainly due to 
several rapidly damped oscillations which, however, are not as well separated 
as in the previous case. 

An integraph? has been constructed which enables the curves for a; and a, 
to be drawn quickly. Fig. 22 is a diagram of this integraph. FF’ is a brass 
block which rolls on two wheels WW’ and carries a pencil at its middle point 
x. The points FF’ on this block are connected to two parallel bars AD and 
A’D’ by means of two Peaucellier linkages which allow F and F’ to move 
along straight lines perpendicular to AD and A’D’ respectively. The bars 
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AD and A’D’ are connected by two links AA’ and BB’ which are supported 
on pivots at S and T fixed to a base. The link AA’ can be turned about S by 
means of a liver SL. This lever also moves a bar ROP which can slide on the 


2 Intergraphs based on the same general principle are described in Les Intergraphs. Ab 
dank-Abakauowiz. Paris, (1886). 
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base parallel to itself. When AA’ is turned the block FF’ also turns so that 
it is always parallel to AA’. When the lever SZ is turned through an angle 6 
the end P of the bar PR moves a distance / tan 6 where / is the distance from 
Sto PR. 

The integraph is moved along parallel to the zero line of the curve of 
which the area is desired and the point P is made to follow the curve by 
moving the lever SL. When the base is moved a distance dx the pencil moves 
along a distance ds such that ds cos 6=dx. Let ds sin @=dz so that dz/dx 
=tan 0. If the first curve is given by y=f(x) then y=/ tan @. Hence y =/dz/dx 
or vdx =/dz so that Is = fydx-s is the distance the pencil moves side ways so 
that the pencil draws a curve giving the area of the first curve. 

In the instrument which was made the lengths 4A’, BB’, FF’ were each 
two inches and / was made one cm. The links were made of stainless steel 
3/8 by 1/16 inch. The wheels were 7/8 inch ball bearings. The block FF’ 
was 3X1X1 inches. 

This instrument enables the curves for a; and a, to be drawn quickly and 
with sufficient accuracy for most purposes. In using it one side of the base is 
kept against a straight edge fastened to the drawing board. 
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The Velocity of Sound in Gases 


N THE classical theory of the propoga- 
| tion of sound in gases it is assumed that 
the gas acts like an elastic medium. The gen- 
eral theory of the propagation of pressure dis- 
turbances (of which sound is one form) 
through such a medium shows that the square 
of the velocity of propagation is equal to the 
rate of change of the pressure of the medium 
with the density (i.e., dp/dp). The formula 
connecting the pressure with the density (or 
volume) of a gas depends upon the supple- 
mentary conditions specifying the state of the 
gas; e.g. whether the process takes place at 
constant temperature (isothermally) or with- 
out the transfer of heat (adiabatically). It de- 
pends also upon the equation of state of the 
gas, of course, but this dependance is not so 
important at temperatures not too near the 
liquefaction point of the gas. By comparison 
with the experimental values for the velocity 
of sound and also by a study of the kinetic 
theory of gases, one concludes that in the 
rapid pulsations occurring in the gas during 


the propagation of sound waves of ordinary ° 


frequencies (40—10,000 cycles/sec.) the gas 
behaves adiabatically as the time interval be- 
tween oscillations is too short for a transfer 
of heat between the various portions of the 
gas. This leads to the formula 


v=[vp/e]? 
where v is the velocity of sound, p the pres- 
sure, p the density, and y is the ratio of the 
specific heat at constant pressure to that at 
constant volume. For simple monatomic gases 
such as He the value of y is about 1.67 while 
for more complicated gases it lies between 1.67 
and 1, being closer to the latter value the more 
complex the structure of the molecule, gen- 
erally speaking. In fact, the law of equiparti- 
tion of energy of Maxwell shows that 

y=1+4+2/s 
where s is the number of independent types 
of motion which the molecule exhibits. While 
this last formula is no longer supposed to be 
altogether correct, if properly interpreted it 
does give one a clue as to the modes of vibra- 
tion of molecule. If one then measures the 


velocity of sound, and supposes its formula as 
given above to be sufficiently substantiated 
theoretically, he can use his measurements for 
a calculation of y, and so can gain some in- 
sight into the complexity of the molecule. 

Recent measurements of this nature have 
been made by H. Kneser' and discussed from 
the point of view of molecular dynamics by 
O. Heil.? Kneser found that while the velocity 
of sound at ordinary frequencies in carbon di- 
oxide gas is about 258.5 meters/second it in- 
creases to about 268.2 meters/second for 
sound frequencies around 10°/sec. This im- 
plies a change in y from about 1.3 to 1.4, and 
a consequent change in s of roughly one unit. 
Presumably then, for very high frequencies 
some mode of vibration of the molecule is not 
excited. 

According to the chemical facts concerning 
CO», the: molecule seems to have a linear 
structure, i.e., the equilibrium position of the 
two O-atoms is on a straight line through the 
C-atom and equidistant from it. The easiest 
mode of vibration to excite in this system, and 
the only one which will be appreciably ex- 
cited at room temperatures, is the oscillation 
of the C-atom normal to the line joining the 
O-atoms, the energy of oscillation of the O- 
atoms themselves being much larger and con- 
sequently harder to excite. It is known, how- 
ever, that it is sometimes quite difficult to 
transfer energy of translational motion into 
energy of vibrational motion in the collision 
between two molecules because of the neces- 
sity for the conservation of momentum as well 
as of energy. In the particular case of CO; it 
appears that only one collision in 10° is effec- 
tive in exciting this vibration of the C-atom. 
Because of this there is not a sufficient num- 
ber of collisions occurring during the time oc- 
cupied by a sound vibration to excite this 
oscillation, and as a consequence the gas acts 
effectively as though its molecules were rigid 
bodies with one moment of inertia (that about 
the common line of the atoms) zero, so that 


1H. Kneser, Phys. Zeits. 32, 179 (1931). 
2 O. Heil, Zeits. f. Physik 74, 31 (1932). 
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the effective number of degrees of freedom of 
the system is 5; i.e., y=1+2/5 =1.4. In meas- 
urements of the specific heat the time re- 
quired for making such a measurement is na- 
turally so long that equilibrium conditions can 
he established within the gas, and this effect 
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does not occur. This intrusion of molecular 
phenomena into a large scale phenomenon 
like the propagation of sound in a gas fur- 
nishes an interesting example of the direction 
of motion of the science of physics. 


eB. L. Fes 





The Pauli Exclusion Principle 


VER since the enunciation of the famous 
KE “exclusion” principle by W. Pauli in 1925 
there has been so much discussion about it,} 
and it has been so fully confirmed from such a 
variety of data that it has become one of the 
most firmly established of physical principles. 
While we cannot hope to bring out and estab- 
lish detailed applications of the principle in 
the following lines, we shall attempt to make 
its meaning a bit clearer to those to whom a 
purely mathematical statement would mean 
little. 

It is important to realize that Pauli ar- 
rived at the statement of his principle from a 
study of experimental data before the de- 
velopment of the quantum mechanics—from 
the data on the character of the orbits occu- 
pied by the electrons in complicated atoms. 
In the atomic theory developed by Bohr the 
attempt was made to assign each electron toa 
certain definite orbit which could be specified 
by quantum numbers. To each orbit there was 
assigned a set of three quantum numbers 2, k, 
and m. From a consideration of the actual way 
in which the data on the spectra of the atoms 
showed these orbits to be occupied, Pauli con- 
cluded that not more than two electrons 
would ever be found to occupy such an orbit. 
The picture of an electron orbit according to 
this theory was an ellipse around which the 
electron travelled rapidly while the nucleus 
was stationed at one focus of the ellipse. The 
integer 2 determined the major axis of the el- 
lipse, & the minor axis, while m determined the 
orientation of the normal to the plane of the 
ellipse in space. 

In order to distinguish between the two 
electrons which could occupy such an orbit 
or “state,” the next step was to introduce a 
fourth quantum number. Thus the state of an 
individual electron is now specified by four 
quantum numbers 2, /, mi, m,; the first three 
have nearly the same significance as in Bohr’s 
theory but the fourth m,, is the so-called spin 
“magnetic” quantum number which assumes 
the two values (+3) just sufficient to dis- 


tinguish the two orbits having the same values 
of n,l, and m). 

It is sometimes said that each electron in 
an atom is specified by certain quantum num- 
bers. This way of speaking is based on the 
older theory in which each individual electron 
was supposed to occupy a given orbit in the 
strict mechanical sense of the word, but it is 
now anachronistic and should not be used 
when there is any danger of being misunder- 
stood. Quantum mechanics as it is under- 
stood at present does not give one the slight- 
est information about what the individual 
electrons may be up to, but instead it at most 
gives one information concerning the possible 
states which the electrons can occupy. It is 
the determination of these allowed states which 
lie at the basis of the theory. These states are 
determined, we suppose, from a “wave equa- 
tion” of some form, either as specified by 
Schrédinger, or perhaps in the (still unset- 
tled) relativistic form of Dirac. 

It is a general rule in quantum mechanics 
(although this point has scarcely been studied 
at all from the mathematical point of view) 
that there are as many quantum numbers for 
the specification of a state of a system as there 
are degrees of freedom (i.e. independent co- 
ordinates) in the system. Since it is found that 
one requires four quantum numbers for the 
complete specification of the individual states 
occupied by an electron in an atom, one would 
suppose that an electron has four degrees of 
freedom, three of which are certainly the three 
degrees of freedom corresponding to the 
Cartesian coordinates of the center. The 
fourth degree of freedom is now interpreted as 
due to the fact that the electron does not act 
altogether like a point charge but has a spin- 
ning motion like a top. In accordance with the 
general principles of quantum mechanics one 
coordinate is required for the specification of 
the direction of the axis of the top in space, 
that is, the space “quantization” of the spin 


1W. Pauli, Zeits. f. Physik 31, 765 (1925). 
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angular momentum vector. This is the rdle of 
the quantum number m, which measures the 
component of the spin angular momentum 
(in units h/27) along some preassigned axis in 
space. In reality a fifth quantum number is 
necessary for the specification of the magni- 
tude of the spinning motion of theelectron, but 
it seems that ordinarily all electrons have the 
same amount of spinning motion around the 
axis of spin, and consequently this fifth degree 
of freedom does not play any important rdle. 
In the relativistic theory of the electron, how- 
ever, the situation is not so simple. 

On the basis of this rather general theory 
it is to be expected, then, that four quantum 
numbers will be needed for the specification 
of the state of an electron even when that elec- 
tron does not form part of an atom; i.e., even 
when it is part of a molecule or of a large 
crystal. This raises the question: will there be 
any analogue of Pauli’s principle applicable to 
electrons which are not part of an atom? 
Clearly there must be, for it can hardly be 
allowed that an atom represents any particu- 
larly peculiar kind of mechanical system. 
What then, is this analogue or generalization 
of the original exclusion principle? 

The answer which theoretical physicists 
have given to this question may be found by 
looking for the essential features of the princi- 
ple as applied to electrons within atoms. This 
generalization is necessarily of a mathematical 
nature in order to be of use in a mathematical 
theory but that does not necessarily mean 
that it is really complex. When we think of 
the states within an atom as being specified in 
some way, through the intermediary use of 
the wave equation, and then of the electrons 
as being assigned to these states, we see that 
Pauli’s principle states in effect: not more 
than one electron is to be allowed to occupy 
a given state (the state being completely speci- 
fied say by the four quantum numbers 1, /, 
m,, m,). This is a general rule which could in 
principle be applied to any problem in which 
the states which electrons can occupy are com- 
pletely specified. As an example we may con- 
sider the problem of the motion of a set of 
electrons within a cubical box which has “re- 
flecting” walls preventing the escape of the 
electrons. Following the original ideas of de 
Broglie, it is convenient to imagine a set of 
waves filling the enclosure. The possible 
“states” for such a system are given by the 
condition that the associated waves shall form 
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a standing wave system. The wave-lengths of 
the waves in each of the three directions nor- 
mal to the sides of the box must accordingly 
be integral sub-multiples of twice the length 
of the box (A=2a/n). Three of the quantum 
numbers defining the states can now be taken 
as the integers 71, m2, m3 specifying the wave 
lengths in these three directions. In order to 
define the state completely we need one more 
quantum number which is to be taken as the 
projection of the spin vector of the electron 
along some arbitrary direction, the quantum 
number m,. This gives us just the model of a 
metal which has been used so successfully by 
Sommerfeld.2 The energy of an electron oc- 
cupying a state specified by the values of the 
quantum numbers 2,=7,'; n2=n2'; n3=n;!; 
m, =m, is found to be 


(h?/8ma?) [(m,")?+ (no')?+ (n3')? | 


and to this approximation does not depend 
upon m,!. 


The generalized Pauli principle would then 
tell us that not more than one electron can 
occupy a state specified by the four quantum 
numbers 11, 2, "3, ms. This restriction would 
be expected to apply regardless of the magni- 
tude of a, that is, regardless of the size of the 
crystal or piece of metal, and is thus not spe- 
cifically a result depending upon the smallness 
of atoms and molecules. A crystal or a piece 
of metal (or even a star) is just as truly a 
single quantum mechanical system as is a 
helium atom. 

While this method of stating Pauli’s prin- 
ciple in terms of states defined by means of 
quantum theory is general, it is by no means 
always convenient, and has been replaced bya 
more succinct and clear statement to the effect 
that thecomplete wave function specifying the 
state of the system (where “the system” now 
means the whole group of electrons and pro- 
tons), must be antisymmetrical in the coordi- 
nates of each pair of electrons— if the coordi- 
nates of any two electrons be interchanged in 
the wave function the latter must be changed 
in sign. In the “coordinates” of an electron 
must be counted a spin variable analogous to 
m, as well as the three Cartesian coordinates 
of the center of the electron. The particular 
advantage of this specification is that one 


2 E. g., A Sommerfeld and N. Frank, Rev. 
Mod. Phys. 3, 1 (1931). 
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need not bother (unless he finds it convenient) 
to specify the actual states which the individ- 
ual electrons are to occupy. 

There is no further theoretical justification 
of the exclusion principle than this. It is not 
a “deriveable” principle—but like wave me- 
chanics itself, draws its justification from its 
agreement with the experimental facts. That 
it has been found to fit in so nicely with wave 
mechanics is, of course, no proof of its validity 
or explanation of its physical content, but 
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nevertheless does give one a more secure feel 
ing concerning it. 

It has been found possible to extend the 
scope of the principle to include the general 
case in which the system contains various 
types of similar particles, such as photons, 
nuclei, etc. For example, it is found that pho- 
tons apparently obey an “inclusion” principle 
in that great numbers of them tend to crowd 
into the same state. 


E. L. Haws 











